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NUMERICAL STUDY OF STEADY AND UNSTEADY 
CANARD- WING-BODY AERODYNAMICS 


Eugene L. Tu 
Ames Research Center 


SUMMARY 

The use of canards in advanced aircraft for control and improved aerodynamic performance is a topic 
of continued interest and research. In addition to providing maneuver control and trim, the influence 
of canards on wing aerodynamics can often result in increased maximum lift and decreased trim drag. 
In many canard-configured aircraft , the main benefits of canards are realized during maneuver or other 
dynamic conditions. Therefore, the detailed study and understanding of canards requires the accurate 
prediction of the non-linear unsteady aerodynamics of such configurations. For close-coupled canards, 
the unsteady aerodynamic performance associated with the canard-wing interaction is of particular 
interest. The presence of a canard in close proximity to the wing results in a highly coupled canard- 
wing aerodynamic flowfield which can include downwash/upwash effects, vortex- vortex interactions and 
vortex-surface interactions. For unsteady conditions, these complexities of the canard-wing flowfield 
are further increased. 

The development and integration of advanced computational technologies provide for the time- 
accurate Navier-Stokes simulations of the steady and unsteady canard-wing-body flowfields. Simula- 
tions are performed for non-linear flight regimes at transonic Mach numbers and for a wide range of 
angles of attack. For the static configurations, the effects of canard positioning and fixed deflection an- 
gles on aerodynamic performance and canard-wing vortex interactions are considered. For non-static 
configurations, the analyses of the canard-wing-body flowfield includes the unsteady aerodynamics 
associated with pitch-up ramp and pitch oscillatory motions of the entire geometry. The unsteady 
flowfield associated with moving canards, which are typically used as primary control surfaces, are 
considered as well. The steady and unsteady effects of the canard on surface pressures, integrated 
forces and moments, and canard-wing vortex interaction are presented in detail, including the effects 
of the canard on the static and dynamic stability characteristics. 

The current study provides an understanding of the steady and unsteady canard-wing-body flow- 
field. Emphasis is placed on the effects of the canard on aerodynamic performance as well as the 
detailed flow physics of the canard-wing flowfield interactions. The computational tools developed to 
accurately predict the time-accurate flowfield of moving canards provides for the capability of coupled 
fluids-controls simulations desired in the detailed design and analysis of advanced aircraft. 




Chapter 1 


Introduction 


1.1 Motivation 

Since the Wright Brothers first introduced the canard configuration almost a century ago, the benefits 
of canards over the now more conventional tailplane configurations have been researched and debated 
extensively. In particular, for the past 30 years, the use of canards in many advanced aircraft for control 
and improved aerodynamic performance has been the topic of continued research. The influence of 
canards on wing aerodynamics can often result in increased maximum lift, and decreased trim drag. 
There are also unique dynamic performance characteristics for canard-configured aircraft coupled with 
the capability of present-day automatic control systems. The reduced or even negative static stability 
of canard configurations can lead to improved aircraft agility and maneuverability. 

There are a variety of aircraft which utilize canards to enhance aerodynamic performance and 
controllability. Aircraft flying today can be found with canards which are fixed or movable, located 
at various longitudinal and vertical positions, and coincident or deflected with respect to the wing. 
For example, the X-31 aircraft has a “no-load” long-coupled canard which deflects between +20 and 
-70 deg for high angle of attack pitch recovery control (ref. 1). Otherwise, during normal flight, the 
X-31’s canard is “free-floating” and deflects to maintain minimal aerodynamic loads. 

Close-coupled canards, by definition, have a more significant effect on the canard-wing aerodynamic 
interaction and, consequently, the aerodynamic performance of the aircraft. The NASA X-29. SAAB 
Viggen and SAAB Gripen are three examples of fully integrated close-coupled canard configurations. 
The X-29 has a forward-swept wing and a movable close-coupled canard which is the primary pitch 
control surface as well as an integral component in the active control system (ref. 2). The SAAB Viggen 
has a close-coupled fixed canard for high performance aerodynamics while its successor, the Gripen, 
utilizes movable close-coupled canards to obtain maximum lift in maneuvering, maximum lift-to-drag 
ratio in cruise, and even nose-down pitching moment during short-field landing roll-out (ref. 3). 

Improved performance of close-coupled canard configurations is usually the direct result of the 
aerodynamic interaction between the canard and wing. However, depending upon geometry and flow 
parameters, this interaction can be either favorable or unfavorable. In many canard-configured aircraft. 
including some of those listed above, the main benefits of the canards are realized during maneuver 
or other dynamic conditions. Therefore, the detailed study of canards as aerodynamic performance or 
primary control surfaces requires the accurate prediction of both the steady and unsteady aerodynamics 
of such configurations. Proper utilization of canards requires the detailed and complete understanding 
of the canard’s steady and unsteady influence on the flow structure about the wing. 
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1.2 Background 


At moderate angles of attack, for aft-swept lifting surfaces with sharp leading edges, the flow can 
separate at the leading edge due to the adverse pressure gradient on the leeward side. Figures 1.1 and 

1.2 show a perspective and crossflow-plane schematic, respectively, of the flow structure over a typical 
canard or wing with leading-edge separation. A free vortex sheet is formed which rolls up over the 
upper surfaces of the canard or wing. If the vortex is strong enough, secondary (S2) and, in some cases, 
tertiary (S3) separations may result. 

The flow structure of highly-swept or delta canard-wing configurations is characterized by a canard- 
downwash which modifies the wing flowfield and an interaction between the canard and wing vortex 
systems. A schematic of a typical canard- wing vortex interaction is given in figure 1 . 3 . In the absence 
of vortex breakdown over the canard, the canard vortex convects downstream and modifies the wing 
flow structure through its induced flow r and interaction with the wing vortex. The downwash of the 
canard modifies the flowfield of the wing within the canard-tip span-line by decreasing the effective 
angle of attack of the wing. Beyond the canard tip, upwash from the canard increases the wing’s 
effective angle of attack. The downwash and upwash effects of the canard significantly influence the 
formation of the wing leading-edge vortex. The canard downwash can weaken or delay the formation 
of the wing vortex, thus affecting its position over the wing surface. 

The effect of changing the position of the canard or deflecting the canard can drastically alter 
this canard-wing aerodynamic interaction. For example, a canard positioned with a vertical offset 
from the w'ing (high- or low-canard) will have varying locations of the canard wake and leading-edge 
vortex relative to the wing, both of which can significantly change the canard-wing flowfield from 
that of the coplanar canard case. For the canard with a positive deflection angle, the stronger canard 
downwash and modified canard trailing edge location will also change the wing flowfield relative to that 
of the coplanar canard case. All of these flowfield changes can potentially affect both the aerodynamic 
performance as well as the stability and trim characteristics of a configuration, and need to be fully 
understood. 

Additional flow features contributing to the complex flow structure of typical canard configurations 
include secondary, trailing-edge, and tip vortices, as well as regions of shock-induced or boundary-layer 
separation. A strong primary vortex often causes the formation of a strong secondary vortex which, 
from the resultant induced velocities, will significantly affect the surface pressures near the canard 
or wing leading edge. Furthermore, trailing-edge and tip vortices can interact with the leading-edge 
vortex as it convects downstream. Boundary-layer separation due to high angles-of-attack or induced 
by a strong recovery shock is also influenced by the presence of these vortices. 

In general, the characteristics of static canard configurations are adequately represented by steady- 
state aerodynamics. At higher angles of attack, some of the conditions which may result in unsteady 
aerodynamics include large regions of separated flow and vortex breakdown. For a configuration un- 
dergoing body motion with a fixed or moving canard, the unsteady effects can be quite significant. In 
particular, the downwash of the canard and the interaction between the canard and wing vortices can 
exhibit highly non-linear unsteady aerodynamic characteristics. 

1.3 Previous Studies 

The use of canards for improved performance has been supported by numerous experimental studies as 
well as some recent computational studies. An early experimental study by Behrbohm (ref. 4 ) indicated 
the potential for improved aerodynamic characteristics of short or close-coupled canard configurations 
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based on the canard- wing interaction. Other experimental studies by Lacey and Chorney (ref. 5), 
McKinney and Dollyhigh (ref. 6), and Dollyhigh (ref. 7) showed the canard’s potential for obtaining 
increased lift and reduced drag. 

Studies by Gloss and McKinney (ref. 8), and Gloss (refs. 9-11) provided further insight into the 
effects of canard geometry and positioning on aerodynamic loading. Experimental studies by Gloss 
et al. (ref. 12). and Dollyhigh (ref. 7) examined the effects of canard size, position and deflection in the 
subsonic and supersonic regimes, respectively. In particular, the studies by Gloss et al. (refs. 10 and 13) 
identified unfavorable characteristics for the low-canard configuration. Substantial loss of lift due to the 
interaction between the low-canard wake and the wing was experimentally measured. However, due to 
the nature of these studies, which only investigated aerodynamic performance parameters, insufficient 
data were obtained to completely describe the canard-wing flowfield. 

Lacey (ref. 14) conducted an extensive four-volume experimental study on the effects of canard 
geometry, position and deflection on aerodynamic loads in the subsonic to supersonic regimes. Other 
studies (refs. 15 and 16) investigated the effects of canard and control surface deflections on more 
realistic canard configurations while references f, 11, and 12 provided detailed studies on generic canard 
configurations with canards at various deflection angles. Recent experimental studies by Howard and 
Kersh (ref. 17), and Ponton et al. (ref. 18) have given extensive information on the flow' structure of 
deflected canard geometries in the low subsonic regime and have shown encouraging results towards 
the optimization of such configurations. 

A series of experimental studies by Gloss and Washburn (refs. 13, 19, and 20) provided detailed 
measurements of surface pressures as well as integrated force quantities on a variety of configurations 
and flow conditions near the transonic regime. Experimental studies by Er-El and Seginer (ref. 21). 
Calarese (ref. 22). and Oelker and Hummel (ref. 23) concentrated on the canard and wing \ortex 
systems and provided details into the mechanisms of their interaction. 

The experimental studies cited above have, in general, concentrated on steady-state aerodynamics. 
Considerably fewer studies have examined the unsteady aerodynamic flowfields, and more specifically, 
the unsteady motion of canard configurations. An earlier study by Bovden (ref. 24) investigated the 
dynamic stability and response characteristics of typical canard configurations and showed potential 
aircraft maneuverability and agility benefits with the use of canards. More recent low-speed water- 
tunnel studies (refs. 25-27) from the Naval Postgraduate School have examined the effects of close- 
coupled canards and canard oscillations on the canard-wing unsteady aerodynamic interaction. A 
low-speed study (ref. 28) of an actual aircraft configuration undergoing canard pitch oscillations has 
also been conducted. 

Computational fluid dynamics (CFD) has become a valuable tool for understanding the complex 
three-dimensional flow physics of canard configurations. In particular, CFD has the potential to greatly 
expand the understanding of the unsteady canard-wing interaction. A number of previous studies based 
on conformal mapping, linear and non-linear vortex lattice methods, the transonic small perturbation 
(TSP) equation, and the Euler equations have been performed (refs. 29-39). However, limited compu- 
tational w r ork has been performed using the time-accurate Reynolds-averaged Navier-Stokes (RANS) 
equations. Although viscous computations are generally not required for certain leading-edge-type sep- 
arations, viscous modeling is essential to capture some of the other significant features of canard-wing 
aerodynamics such as vortex-induced secondary separations and other boundary-layer-type separations. 
With the emergence of faster computers and increased memory capacities, the time-accurate RANS 
equations can now r be utilized. 
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1.4 Objectives 


The primary objective of this study is to develop a complete understanding of the steady and unsteady 
flowfield about canard-wing-body configurations in the high-subsonic and transonic flight regimes. The 
development and enhancement of the computational technology necessary to predict the non-linear 
time-accurate aerodynamics associated with close-coupled canard-wing interactions is an integral and 
necessary portion of this study. In this study, the time-accurate RANS equations are solved for the 
flow about highly-swept, close-coupled canard-wing-body configurations. The major numerical issues 
involved in solving the RANS equations are addressed in detail. 

The accuracy of the computational techniques developed and utilized in this study is demonstrated 
through extensive comparisons with available experimental data. Computational results are compared 
with experimental lift, drag and pitching moment coefficients, as well as surface pressure data. Accuracy 
is further demonstrated through grid refinement studies for the steady-state computations, and through 
both spatial and temporal refinements for the time-accurate computations. Comparisons of the time- 
accurate computational results are also made with the limited unsteady experimental data which are 
available. 

The flow physics associated with canard configurations are studied through extensive computational 
simulations of the steady and unsteady flowfield. The “baseline” geometry with a coplanar canard and 
wing is used to examine the basic steady- flow characteristics of the canard- wing aerodynamic interac- 
tion. These characteristics include the canard-wing leading-edge vortex interaction, canard dowmwash 
effects and wing vortex breakdown. The effects of the canard's vertical position and deflection angle 
on the canard-wing interaction, and subsequent steady-state aerodynamic performance, are investi- 
gated. The time-accurate simulations are utilized to examine the unsteady flowfield characteristics of 
canard-configurations undergoing pitch-up ramp motions, pitch oscillations of the entire configuration, 
and independent canard oscillations. Various pitch rates are used in the study and numerical issues 
affecting the accurate simulation of the unsteady flowfield are addressed. 

In addition to the details of the unsteady flowfield, the influence of fixed and moving canards on 
the dynamic stability characteristics of the configuration are studied. The results of this study not 
only demonstrate the suitability of RANS simulations for aerodynamic analysis and design of canard 
configurations, but also provide for a future direction towards utilizing high-order CFD in stability and 
control simulations. 

1.5 Geometry 

The geometry in this study is illustrated in figure 1.4 and is based on the wind-tunnel models used 
by Gloss et al. (refs. 9, 10, 12, 13, 19. and 20) and Boyden (ref. 24). The current geometry has also 
been tested at low-speeds and high angles of attack by Calarese (ref. 22) in the wind-tunnel facilities 
at Wright Aeronautical Labs. 

The basic geometry consists of a highly-swept wing, close-coupled canard, and a fuselage with a 
constant cross-sectional shape aft of the nose region. The wing and canard are composed of biconvex 
airfoil sections which linearly vary in thickness from 6 percent to 4 percent as defined from root-to-tip. 
The resulting sharp leading edges of both lifting surfaces fix the primary separation locations and are 
conducive to generating vortices at moderate angles of attack. 

For the deflected or moving canard geometries, the canard rotates about the spanwise axis (perpen- 
dicular to the symmetry plane) at the mid-canard-root location as shown in figure 1.4. For the swept 
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canard, static (fixed) or dynamic deflections result in pure rotation at the canard root section and a 
combination of translation and rotation at the canard tip. 

The cross-sectional shape of the body allows for a +10 to -10 deg mid-canard deflection without 
a canard-root /body gap. Canard deflection angles, 5 C , are given relative to the body axis (or wing 
position) with positive 6 C resulting in higher canard angles of attack (a c ). The 6 C = 0.0 deg mid- 
canard deflection represents the coplanar canard case. 



Figure 1.1. Sketch of the leading-edge vortex structure of a canard or wing. 



Figure 1.2. Crossflow plane sketch of the leading-edge vortex. 
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Canard-downwash induced delay 
in wing vortex formation 



Figure 1.3. Schematic of the typical steady-state canard-wing vortex interaction for a con- 
figuration at moderate angles of attack. 


Moment Center 
and Conflg Pitch 



Figure 1.4. Close-coupled canard- wing-body geometry (pitch axis, moment center, and 
canard vertical locations shown). 
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Chapter 2 

Computational Modeling 


The ability to accurately simulate the steady and unsteady flowfield of close-coupled canard config- 
urations is governed by the appropriate computational modeling of the relevant flow physics. Time- 
accurate vortical flows, transonic shocks, and viscous effects including boundary- layer separations are a 
few of the present features which require accurate computational modeling. In addition to the complex 
flowfield. geometric complexities, including moving surfaces, require advanced zonal interfacing and 
dynamic grid generation techniques. 

The computational modeling used in this study is presented in the current Chapter. Geometry 
modeling, grid generation, and other case- dependent issues such as zonal interfacing and boundary 
conditions are discussed in detail. The methods and procedures which comprise the developed com- 
putational technology for accurately simulating canard-wing-body flowfields is described. Details on 
the formulation of the governing equations, numerical algorithm, and turbulence modeling are found 
in the Appendices. 

2.1 Numerical Method 

The numerical method used in this study is based on the solution of the Reynolds- averaged Navier- 
Stokes equations. Due to the significant potential for the presence of transonic and crossflow shocks, 
the strong conservation law form of the N avier— Stokes equations is utilized. As is typical for the 
solution of the governing equations using body-fitted structured grids, the Navier-Stokes equations are 
transformed to and solved in general curvilinear coordinates. 

In general, for high Reynolds number flows, viscous effects are most significant in a thin region near 
solid surfaces. Furthermore, in order to resolve the high gradient flow features normal to the body 
surface, the spatial grid resolution is highest in the body-normal direction. Hence, viscous terms along 
the body are often neglected and the solution of the resulting “thin-layer" Navier-Stokes equations is 
obtained. However, it is noted that for complex configurations, the body normal directions are not 
always confined to one curvilinear coordinate direction. The thin-layer approximation is often made in 
all three directions and. as a result, only the cross-derivative viscous terms are neglected. 

The algorithm employed to solve the Reynolds-averaged Navier-Stokes equations is the Beam- 
Warming alternating-direction- implicit (ADI) algorithm (ref. 40). The form of the Beam- Warming 
algorithm used in the current study is second-order accurate in space and first-order accurate in time. 
Higher-order time accuracy is available with minor modifications, but is generally not required due to 
the small stability-restricted time step sizes. For steady-state applications, the diagonal form (ref. 41) of 
the Beam-Warming algorithm is utilized. The diagonalization of the original block tridiagonal matrices 
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results in scalar pentadiagonal matrices, which have much simpler and more efficient matrix-inversion 
properties with the same second-order spatial accuracy (refs. 42 and 43). The use of spatially- varying 
(or local) time stepping in conjunction with the diagonalization procedure leads to highly efficient 
steady-state computations. Further details on the governing equations and numerical algorithm can 
be found in reference 43 and are also presented in the Appendices. 

2.1.1 Turbulence modeling 

For simulations of complex geometries and flowfields, practical resource considerations dictate the use 
of relatively large space and time scales in the computations. These scales are inadequate for direct 
or large-eddv simulations of turbulent fluid motions. Therefore, the RANS equations are used to 
simulate the flowfield. Using the RANS equations, high frequency small scale motions associated with 
turbulence are not resolved. However, as a result of the RANS "time- aver aging” process, additional 
terms representing effective stresses due to turbulence can be modeled. In this study, the Baldwin- 
Lomax algebraic model (ref. 44) is used to provide this turbulence closure. 

The Baldwin-Lomax model simulates the effects of turbulence by computing an “effective” viscosity 
in the stress terms of the Navier-Stokes equations. This effective viscosity is determined from a two- 
laver algebraic model where the inner layer is a function of density, mixing length, and vorticity 
magnitude. The outer layer is also based on vorticity distribution with additional parameters wffiich 
are functions of wall-normal distance and maximum differences in velocity magnitudes. 

Since the inner algebraic model is a function of vorticity magnitude within the boundary-layer, 
care must be exercised when using the model in vortex-dominated flowfields. In particular, when 
vortices are present near the surface boundary layer, the unmodified Baldwin-Lomax model can yield 
inaccurate turbulent viscosity estimates by scaling with the high vorticity values within the vortex 
core. For the canard and wing with sharp leading edges, these strong vortices are expected to form. 
Therefore, modifications as suggested by Degani and Schiff (ref. 45) are utilized in the current study. 
These modifications provide searching logic within the turbulence model to ensure proper boundary- 
layer vorticity scaling. Further details of the Baldwin-Lomax turbulence model are presented in the 
Appendices. 

2.1.2 Code performance 

Although supercomputing technology is rapidly expanding, particularly in the area of scalable parallel 
computing, code performance continues to be an area of major concern for large-scale complex simula- 
tions. A majority of the CFD community continues to rely on serial supercomputers running efficient 
highly-vectorizable flow solvers. Even though the efficiency “bottleneck” in CFD is often identified 
as geometry definition and grid generation, the ability to compute detailed unsteady flowfields about 
complex geometries is highly dependent on code performance. Memory, CPU cycle requirements and 
convergence or time-step characteristics for a given flow solver are key performance factors. 

The general performance characteristics of the current flow solver including turbulence modeling, 
multi-zoned interfacing (for steady-state simulations), and the application of boundary conditions are 
described below. Approximately 30 words of in-core memory are required for each grid point and the 
typical Cray C90 single-processor performance is measured at 450 MFLOPS and 6 //sec per iteration 
per grid point. It is noted that the performance estimates are dependent on maximum vector lengths 
which is a function of the specific computational-space dimensions. The additional “overhead” for 
dynamic gridding, moving interfaces and more complex boundary conditions, which are required for 
the moving canard cases, increases the cited computational time requirements by up to 100 percent. 
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A significant amount of this increase is due to interface searches and interpolations when mismatched 
moving interfaces are required by geometry complexities. 

For steady-state computations, approximately 4,000 iterations are required to converge from free- 
stream conditions. Convergence is determined by a 3+ order-of-magnitude reduction in L2 Norm 
values as well as constant force, moment and pressure predictions. Angle of attack, Mach number or 
Reynolds number perturbations require an additional 1,500 iterations for re-convergence. The number 
of iterations for convergence increases at higher angles of attack, where vortices are stronger and the 
presence of vortex breakdown is possible. The number of CPU hours required for a typical steady-state 
simulation of the canard-wing-body configuration ranges from 5 to 20 C90 hours depending on grid 
size and convergence characteristics. 

The number of iterations and overall computational time for unsteady cases varies widely. A 
typical unsteady pitch oscillation case requires up to 14,000 iterations per cycle, depending on time- 
step size and grid resolution. For most flow conditions, 2-3 cycles starting from a converged steady-state 
condition are sufficient to obtain a periodic solution. At worst, almost 150 hours of C90 CPU time may 
be required to obtain a periodic condition for a refined-grid simulation of a moving canard configuration. 
Considerably less CPU time is needed for other presented unsteady results using coarser grids, larger 
time steps, and higher reduced frequencies or pitch rates. In fact, key results presented in Chapter 4 
show that higher pitch-oscillation frequencies can be used to significantly reduce the computational 
requirements for dynamic stability predictions. 

2.2 Geometry Modeling 

The fuselage, wing, and canard of the current geometry are analytically defined and, therefore, are 
computationally modeled with a high degree of geometric fidelity. However, in order to simplify the 
modeling, certain aspects of the geometry are approximated. In particular, differences between the 
computational and wind-tunnel models exist in the canard-body junction region, and near the tips of 
both the canard and wing. Other regions which have negligible aerodynamic effects in this study are 
omitted or altered as indicated. 

The geometry used in this study has been tested extensively in a variety of wind tunnels as outlined 
in section 1.5. For these wind-tunnel models, fairings were used to facilitate a vertical-offset canard 
(i.e., non-coplanar) and are present for both coplanar and non-coplanar models. These fairings, which 
account for slight asymmetries in some of the experimental results, are omitted in the current com- 
putational modeling. The vertical tail, which is included in one of the experimental studies (ref. 24), 
is expected to have minimal effects on the static and dynamic longitudinal characteristics and is also 
omitted in the present study. 

The computational modeling for the vertical-offset canard geometry (high- or low-canard) is ac- 
complished by vertically translating the coplanar canard (mid-canard). The surfaces of the canard are 
then extended into the upper or lower portion of the fuselage by linear extrapolations. This extrapo- 
lation results in a thicker canard-root section for the vertical-offset canard. The 6 percent thickness, 
mentioned in section 1.5, applies to the canard span station located at the same span as the wing-root. 

The sting used for wind-tunnel mounting is modeled by extending the body to the downstream 
boundary. In order to compare with the experimental data, integrated force and moment results 
are given for the configuration without the sting and the predicted drag coefficients are corrected to 
free-stream static pressure on the base area of the model. 
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2.3 Static Grid Generation 


Details of the grid generation process for the static geometry simulations are presented in this section. 
A procedure for efficient grid generation of static canard-wing-body geometries has been developed 
and is outlined in figure 2.1. In general for static cases, grid generation occurs only once and can be 
performed as interactive pre-processing or, for sufficiently simple geometries, as an automated initial- 
ization routine within the flow solver code. It is also noted that, as opposed to dynamic geometry cases, 
total computational expenses involved in grid generation of the static geometries are not significant. 

2.3.1 Surface grids 

The baseline surface grids for the static canard-wing-body configurations are generated with algebraic 
techniques. Using the S3D surface geometry and grid generation code (ref. 46) developed at NASA 
Ames Research Center, the canard, wing and body component surface geometries are modeled from 
their original analytical definitions. Figures 2.2 and 2.3 show the baseline and refined surface grids, 
respectively, generated for the wing-body configuration with and without the fixed mid-canard. The 
upper surface of the half-model canard-wing-bodv grid contains 4,625 (baseline) and 8,790 (refined) 
points. To capture complex flow regions, the canard and wing surface grids are clustered near the 
leading edge, trailing edge, and tip. For the wing-body alone geometry, the total number of points on 
the half-body surface is reduced to 2,110 (baseline) and 4,320 (refined) points. 

For the vertical-offset canard geometry, translation, extension and body intersection of the canard 
surface grids are required. These steps are shown schematically in figure 2.4. Starting with the mid- 
canard baseline surface grid, the canard and fore- fuselage are separated from the rest of the surface grid. 
The canard is translated in the vertical direction as required, and then the root section of the canard 
is extended towards the body centerline using linear extrapolation. Finally, surface grid intersection 
techniques employed in the S3D (ref. 46) program are used to produce a single, continuous, offset 
canard-wing-bodv grid. An example of the surface grid for the high-canard case is shown in figure 2.5. 

Deflecting the canard is performed with an algebraic redistribution of the undeflected canard surface 
grid combined with cubic spline interpolation. As with the vertical-offset canard geometries, only the 
canard portion (canard plus forebody) of the surface grid is modified for the deflected canard cases. 
The generation of the deflected canard surface grid, shown schematically in figure 2.6, is outlined as 
follows: 

1. Starting with the undeflected canard grid, the canard points are rotated about the specified axis; 

2. The upper/lower body dividing line forward of the canard is assumed to be parabolic; 

3. The axial distribution of the body cross-sections is adjusted according to the rotated canard: 

and 

4. The points along each cross section are redistributed using cubic spline interpolation based on 
the new arc length of each section. 

It is emphasized that for all deflected canard cases, the canard is actually rotated (not sheared) 
about the mid-canard-root spanwise axis to maintain high geometric fidelity. An example of the surface 
grid for the deflected canard case ( 6 C = 10.0 deg) is shown in figure 2.7. 

2.3.2 Flowfield grids 

The static flowfield grids are generated by solving the elliptic grid generation equations with orthonor- 
mal and clustering control functions near the zonal boundaries. Such techniques are readily available 
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in grid generation programs such as 3DGRAPE (ref. 47) and GRIDGEN (ref. 48). both of which are 
utilized for the current canard-wing-body geometry. 

The baseline grid for the undeflected mid-canard configuration is shown in figure 2.8 and is generated 
as a single block, H-0 topology grid with 152 axial. 32 radial, and 97 circumferential points. The H-0 
topology is ideal for this type of configuration because of the natural clustering of grid points in certain 
regions of the grid. When clustering is applied to obtain a viscous grid at the canard, wing, and body 
surfaces, grid clustering also results along the complex flow regions of the canard and wing wakes, and 
in the singularity regions along the nose radial axis and upstream of the canard. In order to resolve 
the viscous boundary layer near the canard, wing and body surfaces, an average value of y~ < 3 for 
the first point away from the wall was achieved by clustering the flowfield grid to the surface (average 
first grid spacing is approximately 10 _6 c). where y + is defined as (ref. 49) 



In equation 2.3-1. y is defined as the wall-normal coordinate direction. The stretching of the grid in the 
wall- normal direction is also minimized (typically, stretching factor < 1.2) to reduce numerical errors. 

For the refined flowfield grid, the total number of points is increased from approximately 470.000 
to over 1.7 million. As shown in the expanded view of figure 2.9, refinement of the grid is emphasized 
in several key regions of the flowfield. In addition to an improved surface grid distribution, refinement 
is achieved in the interaction region between the canard and wing. The viscous regions near all of the 
surfaces are also refined by reducing the stretching factor while maintaining the original y + criteria. 
By reducing the amount of grid stretching, the resolution of any vortices near the canard and wing- 
surfaces will be improved. 

Earlier figures 2.5 and 2.7 showed that the generation of vertical-offset, or deflected canard surface 
grids produces a mismatch on the body at the canard-wing junction. Furthermore, a geometric mis- 
match is observed between the trailing-edge root of the close-coupled canard and the leading-edge root 
of the wing. Consequently, to maintain the desired H-0 topology grids, a mismatch is produced in the 
zonal interface separating the canard and wing portions of the flowfield grid. 

The flowfield grids of the different canard position and deflection cases are accommodated by 
generating the canard and wing portions of the flowfield grid individually, with a mutual overlapping 
interface analytically defined. Then, as required, the flowfield grids for the vertical offset are re- 
generated using the elliptic solvers. 

For the deflected canard cases, rather than using interactive elliptic solvers, the flowfield grids are 
re-generated using an algebraic redistribution of the original elliptically-generated coplanar canard- 
wing-body grid. Figure 2.10 illustrates the simple technique used in the algebraic redistribution. The 
redistribution is performed for each line in the £ computational direction (from the surface to the outer 
boundary). Based on the new position of a given surface point, each off-body point is repositioned using 
cubic spline interpolation scaled with arc length from the outer boundary. Although this procedure 
does not guarantee grid quality for large surface grid motions, it produces good quality grids and 
is extremely fast and robust for small surface perturbations. As such, the highly efficient algebraic 
redistribution procedure is well suited for the unsteady moving canard cases described in the next 
section and Chapter 5. 

By using these techniques for both the surface and flowfield grids, the grids for any canard vertical 
offset and deflection angle (6 C ) can be efficiently obtained. The overall H-0 topology grid, with the 
mismatched interface, for a typical vertical-offset and deflected canard case is shown in figures 2.11 and 
2.12. respectively. An expanded crossflow plane view of the mismatched interface is also shown with 
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the solid grid indicating the downstream boundary of the canard portion grid. The zonal interfacing 
techniques used to handle the mismatched interfaces are described in section 2.6. 

The number of points for the resulting flowfield grids of the coplanar, vertical-offset and deflected 
canard configurations are summarized in table 2.1. Where needed, up to six zones stacked in the 
streamwise direction are used to minimize core memory requirements. Except for the canard-wing 
interface of the vertical-offset and deflected canard cases, one-to-one grid matching is used on all 
interfaces. Since the current computations are performed in the high-subsonic and transonic regimes, 
all of the flowfield grids are extended upstream and downstream by approximately eight wing root-chord 
lengths, and in the radial direction by six wing span lengths. 


2.4 Dynamic Grid Generation 

The generation of surface and flowfield grids for moving geometries can be a complex and time- 
consuming process. In particular, for motion which is not pre-determined (e.g., dynamic maneuvers, 
free aeroelastic responses, active control surface motions), grid generation is required to be an integral 
part of the time-iterative flow solution. For prescribed motions such as forced oscillation, the process 
can be simplified. With sufficient computational storage resources, the gridding of predescribed mo- 
tions can be obtained by generating the grids as a pre-processing step to the flow solver. However, 
this storage requirement becomes prohibitively large for even simple cases requiring viscous analyses. 
Therefore, the capability to integrate efficient surface and flowfield grid generation within a flow solver 
is highly desirable and is often needed. 

Figure 2.13 illustrates the integration of the numerical flow solver and dynamic grid generation. 
As each step of the time-integration is completed, the new surface geometry position is determined. 
For the cases presented in this study, the configuration and canard motions are based on prescribed 
ramp motions or forced oscillations. However, it is noted that for dynamic maneuver simulations, the 
equations governing the dynamics of the configuration can be easily integrated into the flow solver/grid 
generation process described in this study. 

Based on the new geometry position, surface grids are generated accordingly. For the moving 
canard configuration of the current study, the steps outlined in section 2.3.1 for deflected canards are 
utilized. The flowfield grid is then generated based on the algebraic redistribution concept described in 
section 2.3.2 and illustrated in figure 2,10. By using these techniques, dynamic grid generation of the 
moving canard flowfield produces less than 10 percent in grid generation overhead with no additional 
computational storage requirements. 


2.5 Boundary Conditions 

The solution of the Navier-Stokes equations requires boundary conditions at each of the block bound- 
aries. Although the finite difference scheme outlined in section 2.1 is implicit, the boundary conditions 
are applied explicitly in the current study. Numerous past computational studies have demonstrated 
the successful coupling of explicit boundary conditions with implicit flow solvers. 

The boundary conditions which are required for external flow about canard-wing-bodv configura- 
tions can generally be categorized into two main types: 1) physical boundary conditions such as solid 
surfaces or outer boundaries of the solution flowfield, and 2) zonal boundaries or interfaces between 
adjacent (or the same) computational blocks. The first type of boundary condition is described below 
in this section. The implementation of the zonal boundary or interfacing is presented in section 2.6. 
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Physical boundary conditions are applied to the boundaries of the entire computational domain. 
For the current canard- wing-body study, these boundaries include the solid surfaces of the configuration 
(canard, wing, and body), inflow boundary, outer computational surfaces, and the outflow boundary. 
The solid surface boundary conditions are generally referred to as “no-slip” conditions. The components 
of flow velocity in physical Cartesian space are explicitly set to the physical velocity of the solid surface. 
The resulting types of solid surface boundary conditions for velocity are shown figure 2.14 and are 
outlined as follow's : 

Type 1) For a static configuration, the physical and contravariant velocity components at solid 
surfaces are set to zero; 

Type 2) For a configuration with a fixed canard in rigid body motion, the surface grid remains 
fixed with respect to the physical surface and the physical velocity is set to the surface velocity (the 
contravariant velocity remains zero); 

Type 3) For the moving canard configuration, the surface grid of the body (fuselage) wall be moving 
relative to the physical surface. In this case, the surface velocity is computed separately from the grid 
velocity and is used to set the physical velocity components. If the entire configuration and the canard 
surface are in independent motion, both the physical and contravariant components of velocity will be 
non-zero. The additional computational step required for Type 3 boundary conditions is the separate 
evaluation of the surface velocity at a physical point from the evaluation of the total grid velocity. 

Surface pressures and densities are determined with the aid of the normal momentum equation, 
perfect gas assumption, and an adiabatic wall condition. For flat surfaces, the normal momentum 
equation yields the familiar dp/dn = 0. For the more general case, terms due to wall curvature and 
acceleration of the moving surface yield 


dp 

dn 


p u 

~r 


+ pdf, ■ n 


(2.5-1) 


where R is the local wall curvature and (a b - n) represents the acceleration of the physical point on the 
surface in the outward normal direction. The first and second RHS terms in equation 2.5-1 represent 
effects due to wall curvature and wall acceleration, respectively. 

However, further investigation of equation 2.5-1 reveals that differences in pressure due to both 
wall curvature and surface acceleration effects scale directly with the first spacing of the wall-normal 
grid. For viscous grids, the first spacing is small enough such that dp ~ 0 is a good assumption and 
zeroth-order extrapolation of pressure is reasonable. Density at the surface is determined by assuming 
an adiabatic wall condition which also yields a zeroth-order extrapolation of both temperature and 
density. 


2.6 Zonal Interfacing 

The close-coupled nature of the canard and wing surfaces can produce significant geometric and grid 
topological discontinuities for deflected or vertical-offset canard configurations. As described in sec- 
tion 2.3.2 and in figures 2.11 and 2.12, a mismatched interface in the flowfield grid can result between 
the canard and wing zones. 

For the mismatched interface in this study, bilinear interpolation as formulated in references 50 
and 51 is used to transfer flow quantities from one zone to the other. For general cases where the 
mismatched interface is not necessarily composed of coincident surfaces, trilinear interpolation may 
be required. Since this interpolation is non-conservative, accuracy has been determined by comparing 
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the solutions from a single-block (conservative) case with the multi-block mismatched-interface (non- 
conservative) case for the undeflected mid-canard configuration. Differences in surface and flowfield 
pressure contours between the two test cases were determined to be negligibly small. 

As illustrated in figure 2.15. the mismatched interface between the canard and wing is overlapped 
by one streamwise grid cell and actually consists of two sets of “base” and “target” surfaces. The 
target surface represents the outer boundary of the current zone being solved, while the base surface is 
a subset of an adjacent zone and is used to provide flow information from that zone. The interpolation 
is performed by “locating” each point of the target surface within the set of base surface points. The 
search required to identify the four base points which comprise the “cell” containing the target point 
can be quite time consuming. After the target point is located, weighted coefficients (representing 
distances) are computed for the nearest three base points. These coefficients are then easily used to 
bilinearly interpolate the appropriate flow quantities to the target surface. 

For static canard configurations, the interpolation described above is required only once and the 
weighted coefficients are stored for use throughout the solution process. However, for the moving canard 
cases, this mismatched interface will change for each time iteration and the interpolation is required 
for each iteration. Figure 2.13 provides a flow-chart diagram illustrating the integration of the zonal 
interface procedure within the solution process. After the new surface and flowfield grids are generated 
for the current zone and time iteration, the interpolation for the mismatched interface is performed 
and the appropriate flow quantities are transferred to adjacent zones. Since this process is required for 
each time iteration, local search techniques which can significantly improve the interpolation efficiency 
are employed. 

Figure 2.16 illustrates a search procedure which, given a favorable starting position (cell #1), can 
significantly reduce the time required to perform the interpolation. For the current study with moving 
grids, a starting position based on the given target point position of the previous iteration is used. 
Since time steps and the corresponding grid movements are relatively small, a successful search and 
interpolation is usually accomplished within 10 search steps. The procedure is then repeated for each 
target point in the target surface. Using this technique, the time required to interpolate the interface is 
reduced by up to two orders of magnitude. For the moving canard cases, interface interpolation results 
in a 25-50 percent approximate increase in overall computer time requirements. 


2.7 Dynamic Stability Analysis 


In order to understand better the unsteady characteristics of the canard-wing-body configurations, 
dynamic stability analyses are performed on the time-accurate computational simulations. Both rigid 
pitch-oscillation of the entire configuration and independent canard oscillations (moving canard cases) 
are considered in the analyses. The effect of the presence of the canard on dynamic stability parameters, 
including pitch damping and oscillatory stability, are considered. 

Using linear theory, the time-dependent moment coefficient due to rigid pitch oscillation of the 
entire configuration can be expressed as 


Cm — -}- 


ac 

2UZ 


C m . + 


_qc_ 
2U « 


Crrin T 


qc~ 

4U 2 

oo 


C mA + 


... higher order terms (h.o.t.) (2.7 — 1) 

where C mQ , C md , C m<7 , and C m - are the stability derivatives (dC m /da) o, (dC m /d(^r^))o 

■ -2 

{dCm / o and (dC m /d(-^r))o. respectively; q is the angular velocity of the configuration about 
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the pitch axis; and C mo = 0 (not included in eq. 2.7-1) for the symmetric configuration. The notation 
()o indicates that partial derivatives are evaluated assuming no disturbance from the other terms. 

Assuming a sinusoidal motion of a = a m + Asinut, where A is the amplitude of the oscillation, 
equation 2.7-1 can be expressed as 

Cm z ot m Cm a T ( Cm<x k C m ^)AsinLutA' 


( Cma T C mq ) Akcosujt + ...h.o.t. (2.7 2) 

The terms ( C ma - k 2 C m and (C ma + C m<t ) in equation 2.7-2 represent the oscillatory stability and 
damping-in-pitch parameters, respectively. 

For the oscillating canard cases, the sinusoidal canard deflection as a function of time is given as 
<5 C = 6 cm + Acsirvju'ct , where A c is now the amplitude of the canard oscillation. In this case, the 
time-dependent moment coefficient due to canard oscillations is written as 


Cm — OcmC i 


cm^mj 


+ 


2 U 0 


-C r , 


+ 


qcc 

2f7oc 


C„ 


• -9 

<? C C~ 

4C/2 


Cm,**' T ...h.o.t. 


(2.7-3) 


where C Tn . . C m i C ma > an d C m „ are the stability derivatives (dC m /d6 c )o. ( dC m jd[ ) ) o 
{dCm!d {^)) o and (dC m /d(ffjr))o- respectively; q c is the angular velocity of the canard about its 
pitch axis. 

Again, the resulting dynamic stability parameters are derived from 


Cm — fi Co C m6c + ( C m6c k c Crn qc )A c sinuJ c t-\- 

( Cm 6c + C mqc )A c k c cosuj c t + ...h.o.t. (2.7 - 4) 

where the stability terms are now expressed as (C mSc — k 2 C mqc ) and (C m ^ + C rnqc ) representing the 
oscillatory stability and damping-in-pitch parameters, respectively, for canard oscillations. 

Similar analyses of the normal force coefficient given by 


C.y = QCy Q + 


ac 

2U X 


C.\ r ; + 


QC 

2 U a 


■-2 

r i 9C 

jV ’ + 4 Ul 


C Nq ■+■ ...h.o.t. 


(2.7-5) 


or 
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(2.7-6) 


will yield the pitch displacement and pitch rate terms given by (Cv Q - k 2 C n-) and (C.y Q + Cjv,) from 
eq. 2.7-5; or iC — kc 2 C ,\ qc ) and ( Cn ^ + C,v ?c ) from eq. 2.7-6. 

The dynamic stability parameters for the canard configurations are predicted from the Fourier 
analysis of the computed time histories of normal force and pitching moment coefficients. For example, 
the time-history of computed pitching moments can be expressed as 


Cm(t) = 


do 


4- aisiruot + bicosiot + ...h.o.t. 


(2.7-7) 


where the Fourier coefficients are given by 
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and 


2 f T o/ 2 

b n = — / C m (t)sin(nu>t)dt ; n — 1,2, ... 
io J —T 0 /2 

In most cases of forced sinusoidal oscillations of the canard or the entire configuration, the normal 
force and pitching moment responses are adequately represented by the first three Fourier coefficient 
terms (ao, ai, and &i). This assumption is easily verified by direct comparisons of the original predicted 
responses and the Fourier analysis using only these terms. 

A direct comparison of equation 2.7-7 with equations 2.7-2 or 2.7-4 yields the result 

{Cm a — k 2 Cm q ) or ( Cm 6c — k?C mqc ) = (“ 7 t ) or ( * , ) 

for pitch damping, and 

(Cm a + Cm q ) or (C m , e + C m J = (^) or (^) 

for pitch oscillatory-stability parameters. Similar analyses of the normal force responses will produce 
predicted values for the pitch displacement and pitch rate parameters. 


Table 2.1. Surface and flowfield grid data for coplanar, vertical-offset and deflected canard configurations 


Configuration 

Surface 

points 

(half-body) 

Flowfield grid 
dimension 
( j , k, 1- nor m al ) ( zone ) 

Total 
number of 
grid points 

Baseline canard-off 

4,220 

99 x 79 x 32 

250.000 

Refined canard-off 

8,640 

136 x 119 x 58 

940.000 

Baseline mid-canard 

9,250 

152 x 97 x 32 

470.000 

Refined mid-canard 

17,580 

250 x 119 x 58 

1.700,000 

Low- or high-canard 

11,000 

89 x 128 x 40 (canard) 
70 x 97 x 40 (wing) 

725.000 

Deflect canard 

9,250 

84 x 97 x 40 (canard) 
70 x 97 x 40 (wing) 

600.000 

Unsteady fixed canard 

9,250 

152 x 97 x 32 

470.000 

ramp or pitch (baseline) 

Unsteady fixed canard 

17,580 

250 x 119 x 58 

1,700.000 

ramp or pitch (refined) 

Unsteady moving canard 

9,250 

84 x 97 x 40 (canard) 
70 x 97 x 40 (wing) 

600.000 
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Figure 2.1. Schematic diagram of the grid generation procedure for static canard configu- 
rations with various canard vertical positions and deflections. 
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Figure 2.2. Baseline surface grid for the wing-body configuration with and without an 
undeflected mid-canard. 




Figure 2.3. Refined surface grid for the wing-body configuration with and without an 
undeflected mid-canard. 
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Figure 2.4. Schematic of the procedure for generating a vertical-offset canard configuration 
from a baseline undeflected mid-canard configuration. 



Figure 2.5. Perspective view of the surface grid for wing-body configuration with a high 
canard. 
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Figure 2.6. Schematic of the procedure for generating a deflected canard configuration from 
a baseline undeflected mid-canard configuration. 


MISMATCHED SURFACE GRID 



Figure 2.7. Perspective view of the surface grid for wing- body configuration with a deflected 
canard (£ c = 10°). 
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Figure 2.8. Flowfield grid topology and expanded near-body grid for the undeflected mid- 
canard configuration. 
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BASELINE GRID 
152 X 32 X 97 


REFINED GRID 
250 X 58 X 119 


Figure 2.9. Comparison of the baseline and refined near- body ilowfield grids for the unde- 
flected mid-canard configuration. 
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(S)k 

(S)k=1 


Figure 2.10. Algebraic redistribution of the flowfield grid. 
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Figure 2.11. Flowfield grid at the mismatched interface surface for the high-canard config- 
uration. 



Figure 2.12. Flowfield grid at the mismatched interface surface for the deflected canard 
configuration. 
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Figure 2.13. Schematic illustrating the integration of dynamic grid generation and zonal 
interface computation into the time-iterative solution process. 
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TYPE 1 

NON-MOVING WALL 
NON-MOVING GRID 
PHYSICAL VELOCITY = 0 
CONTRAVARIANT VELOCITY = 0 



a) Fixed canard steady-state case 


TYPE 2 

MOVING WALL 
GRID MOVING WITH WALL 
PHYSICAL VELOCITY * 0 
CONTRAVARIANT VELOCITY = 0 



TYPE 3 

STATIC/MOVING WALL 
MOVING GRID INDEP. OF WALL 
PHYSICAL VELOCITY = 0 or * 0 
CONTRAVARIANT VELOCITY * 0 





c) Moving canard unsteady case 


Figure 2.14. Types of velocity boundary conditions at the surfaces of the canard, wing and 
body. 
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Interface 

Surfaces 


Figure 2.15. Schematic illustrating the interfacing of two adjacent zones. 



Figure 2.16. Example of a local search procedure for improved efficiency in dynamic zonal 
interface computations. 
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Chapter 3 

Effect of Canard on Steady-State 
Aerodynamics 


The influence of canards on the steady-state aerodynamics of the wing-body geometry is studied 
through a series of computational simulations. Computations are conducted for the configuration with 
and without the canard, with the canard at various vertical positions, and with the canard at specified 
deflection angles. The accuracy of the steady-state computations is assessed with experimental com- 
parisons and a grid refinement study. The effect of the canard is characterized through analyses of the 
aerodynamic performance (e.g., lift, drag, and pitching moments), and with the detailed investigation 
of the canard-wing flow structure. 


3.1 Experimental Validation: Coplanar Canard 

To assess the accuracy of the computational modeling given in Chapter 2 for steady-state predictions, 
computations are conducted on the wing-body geometry with and without a coplanar canard. All 
computational results in this section are for a transonic Mach number of 0.90, a Reynolds 

number based on mean wing aerodynamic chord {Re~c) of 1.52 million, and nominal angles of attack 
(a) ranging from 0 to 12 deg. Based on transition strips used in the experimental studies, boundary- 
layer transition for the computations is assumed to occur at the canard and wing leading edges and 
at the nose of the body. Comparisons between the computed results and the force balance and wing 
surface pressure measurements are made for both the baseline and refined grids. 

3.1.1 Baseline grid 

The first set of results is for computations performed on the baseline grids which were generated from 
the surface grids of figure 2.2. A comparison of computed wing surface pressure coefficients ( C p ) with 
experimental data (ref. 13) at a « 4 deg is illustrated in figure 3.1. The precise angles of attack are 
reported in reference 13 and are shown in the figure. It is noted that no angle-of-attack or Mach number 
corrections to calibrate the code were made for any computations in this study. Comparisons with the 
experimental data are given at wing semispan-stations of 25, 35, 45, 65. and 85 percent as measured 
from the symmetry plane. The first three span-stations on the wing are within the canard-span region. 

For the canard-off case, the leading-edge vortex, as indicated by the suction peak in the C p distri- 
bution. moves aft with increasing wing span. At the relatively low incidence of approximately 4 deg, 
the major canard influence on wing surface pressures is the canard downwash effect. Within the span 
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of the canard, a significant difference in wing surface pressure between the canard-on and canard-off 
cases is observed. The effective local angle of attack (a e ff) of the wing is reduced and the formation 
of the wing leading-edge vortex is inhibited. At wing span-stations beyond the canard tip, the wing is 
no longer influenced by the canard downwash and the leading-edge vortex forms. In fact, outboard of 
the canard tip the canard upwash increases the Oc e ff of the wing. This normally would be expected to 
increase the lift of the wing-tip region. However, the upwash effect is minimized due to the high sweep 
of the wing, which results in a relatively large distance between the canard-tip and wing-tip regions. 

The discrepancies between the computed and measured data towards the wing tip are due to both 
computational modeling of the tip geometry and grid resolution. A previous study by Srinivasan et al. 
(ref. 52) found that, especially for separated and vortex-dominated flows, the modeling of the tip 
geometry has a considerable effect on the accuracy of the numerical results near the tip. Since detailed 
geometric data were unavailable for the canard and wing tips of the wind-tunnel model, the tips are 
computationally modeled as rounded. 

In addition to the tip modeling, grid resolution has a significant effect on the accuracy of the 
computations. Difficulty in capturing the leading-edge vortex increases towards the canard and wing 
tips. As the vortex is convected downstream and away from the leading edge, computational accuracy 
becomes increasingly sensitive to the degrading grid resolution. This effect becomes more acute at 
higher angles of attack where the vortex trajectory is further above the wing surface. An effort to 
improve the computational accuracy by increasing the grid resolution is performed and presented later 
in this study. 

Although detailed comparisons of flow quantities such as surface pressure are better measures of 
the computational accuracy, integrated quantities such as lift, drag, and moments are often used to 
assess the overall aerodynamic performance characteristics of a given configuration. Figure 3.2 presents 
the comparison of canard-on and canard-off integrated force and moment quantities for the coplanar 
canard- wing-bodv configuration. The canard-on computations are also compared with the appropriate 
experimental data (ref. 10). Figure 3.2 shows that the nonlinearity of the lift coefficient curve for the 
canard-on case is captured well by the computations. For angles of attack less than 6 deg, the computed 
lift coefficients for the canard-off and canard-on cases are comparable. At these low angles of attack, 
the reduction in wing lift due to the canard downwash is balanced by the additional canard lift. As 
the angle of attack is increased, the canard-on lift curve exhibits significant nonlinearity. 

The drag polar in figure 3.2 indicates comparable levels of drag coefficient for both the canard- 
on and canard-off cases. The "cross-over” of the two drag curves shows the potential of the canard 
configuration for reduced drag at a given lift. At low angles of attack, the canard-on drag is higher due 
to the additional viscous drag of the canard surface area. At higher angles of attack, the computed 
additional lift due to the canard and the canard's influence on the wing aerodynamics results in a 
higher lift-to-drag (L/D) ratio for the canard-on case. 

Due to the relative location of the canard and the moment-center (fig. 1.4), the pitching moment 
curve in figure 3.2 illustrates the typical nose-up pitching moment which is characteristic of many 
canard configurations. At higher angles of attack, the computations overpredict the pitch-up moment 
of the canard-on case. This overprediction is further studied by examining the lift and moment curves 
of the canard and wing regions separately. 

The lift and pitching moment coefficient curves for different regions of the canard- wing-bodv config- 
uration are given in figure 3.3. The forward and aft regions including the body (henceforth designated 
as the canard and wing regions, respectively) are chosen to correspond with regions measured in the ex- 
periment (ref. 10). The significant decrease in the wing-region lift for the canard-on case in comparison 
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with the canard-off case is due to the canard downwash. In figure 3.3. the lift coefficient of the com- 
plete canard-on configuration is being supplemented by the canard region lift. Favorable experimental 
comparisons of the canard and wing region lift curves confirm that the computed lift is accurately 
distributed between these two regions. The pitching moment curves show that the canard's influence 
on the wing region pitching moment is nominal. The nose-up pitching moment of the canard-on config- 
uration is almost entirely due to the canard region. Again, good comparisons with experimental data 
indicate accurate distribution of pitching moment in the computed results. Figure 3.3 also shows that 
the overpredicted pitch-up moment for the complete configuration, noted in figure 3.2, is approximately 
evenly distributed between the canard and wing regions. 


3.1.2 Refined grid 

At higher angles of attack, the canard and wing leading-edge vortex trajectories move above and away 
from the highly-clustered grid region of the canard and wing surfaces. In addition, separated regions on 
the wing become more prominent, recovery shock strength increases and the potential for leading-edge 
vortex burst exists. Accurate computational modeling of these and other relevant flow features requires 
improved surface and flowfield grid resolution. 

To improve the accuracy of the computed results in this study, refined surface and flowfield grids 
were generated for the coplanar canard-on and canard-off cases. The refined surface grids were shown 
earlier in figure 2.3. Refinement of the wing and canard surface grids, as well as the body surface grid 
near the canard-wing junction, was emphasized. For the current topology, the body surface grid at 
the junction directly determines the grid resolution between the canard and the wing, and therefore 
affects the computational accuracy of capturing the canard's influence on the wing. The resulting 
canard-off and canard-on flowfield grids generated from the refined surface grids contain over 900.000 
and 1.7 million points, respectively. 

Since the wing, in the absence of the canard, exhibits a higher o e // and a more pronounced 
leading-edge vortex, the canard-off configuration w r as chosen to verify the anticipated improvement in 
computational accuracy with the refined grid. Figures 3.4 and 3.5 illustrate computed surface pressure 
coefficients for the baseline and refined grids of the canard-off configuration at a = 8.21 deg and a = 
12.38 deg. respectively. In both figures, the comparisons with the experimental data show considerable 
improvement for the refined grid solution over the baseline grid solution. Figure 3.6 illustrates the 
computed force quantities using the refined canard-off grid. Note that since the computations were 
performed at three angles of attack only, the refined grid results in figure 3.6 are indicated as solid 
square symbols. Improvements in lift, drag and pitching moment comparisons with experimental data 
are also observed. 

A comparison between canard-on and canard-off wing surface pressure coefficients at a « 12 deg is 
shown in figure 3.7. Good agreement between the computed and experimental data (ref. 13) shows the 
suitability of both the canard-on and canard-off refined grids for angles of attack up to at least 12 deg. 
A strong secondary vortex is evident from the sustained low pressure region between the suction peak 
of the primary vortex and the wing leading edge. In comparing the results of figure 3.7 with figure 3.1, 
one observes that for the canard-on case at higher angles of attack, the wing leading-edge vortex is no 
longer inhibited at the inboard wing span-stations. The vortex forms at the wing apex and remains 
near the leading-edge up to the wing location corresponding to the canard-tip span-line. 
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3.2 Aerodynamic Performance 


The highly interactive flowfield of a close-coupled canard results in aerodynamic performance which is 
extremely sensitive to changes in canard geometry, position and deflection angle. These changes, in 
turn, have pronounced effects on the lift, drag and pitching moment performance of the wing and the 
entire configuration. In order to determine the potential aerodynamic performance benefits of utilizing 
canards, the effects of canard position and deflection on aerodynamic performance are investigated in 
this section. 

Comparisons between the computed results and the available experimental data (refs. 9. 10. 12. 13, 
19, and 20) are made to validate the current method for accurately predicting the stead}' flowfield of the 
vertical-offset and deflected canard cases. Computations for these cases are at transonic Mach numbers 
(Mqo) of 0.85 and 0.90, angles of attack (a) ranging from -4 to +12 deg, canard vertical positions of 
z/c = -0.185, 0.0 and +0.185 (henceforth referred to as low-, mid- and high-canard, respectively), and 
canard deflection angles (<5 C ) for the mid-canard cases from -10 to +10 deg. To match experimental 
conditions, at = 0.85 and 0.90. the Reynolds number based on the mean aerodynamic chord of the 
wing (Rec) is 2.82 million and 1.52 million, respectively. As with the earlier coplanar (or undeflected 
mid-canard) cases, transition strips were used in the experiments at the body nose and leading edges 
of the canard and wing. Therefore, the computations are performed assuming fully turbulent flow. 

3.2.1 Canard vertical position 

Comparisons of the lift, drag and pitching moment curves with experimental data (ref. 10) are given 
in figure 3.8 for the low-, mid-, and high-canard cases. Figure 3.8(a) shows the nonlinearities in the 
high- and mid-canard lift curves. At a > 4 deg, significant differences are noted in the lift curves 
for the three canard positions. Although the absolute values of lift are slightly underpredicted by the 
computations, the relative differences are captured. Further investigation of the incremental changes 
due to canard vertical position is presented later in this section. 

Figures 3.8(b) and (c) show computed drag and pitching moment coefficient comparisons with 
experimental measurements. At the higher angles of attack, the low-canard configuration exhibits lower 
drag at a given a. Figure 3.8(c) shows an overprediction of pitching moment by the computations at 
the higher angles of attack. This overprediction is of the same order-of-magnitude as shown earlier and 
discussed in figure 3.2. However, as in the lift and drag curves of figures 3.8(a) and (b), the overall 
pitching moment comparisons are favorable. 

Figures 3.9(a) and (b) illustrate lift and moment curves, respectively, for the two component regions 
of the high-canard configuration. As with the coplanar canard studies, the canard region consists of 
the canard and the body forward of the wing leading-edge root location (fore-body). The wing region 
consists of the wing and the remaining aft-body (not including the sting). 

The wing portion lift curve slope of figure 3.9(a) is influenced by the effect of the canard downwash 
on reducing wing lift. The total lift is augmented by the canard portion lift. Figure 3.9(a) also shows 
that the slight underprediction in total lift is primarily due to the computed wing portion lift. The 
pitching moment curves in figure 3.9(b) show that the wing, in the presence of the canard, exhibits 
statically stable characteristics while the canard pitching moment causes the total configuration to 
be approximately neutrally stable. Such characteristics are often typical of canard-configured aircraft 
and were also evident in the mid-canard configurations shown earlier. Favorable comparisons in fig- 
ures 3.9(a) and (b) between computations and experimental data indicate that the distribution of lift 
and moment between the component regions is captured well by the computations. 
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In order to predict the potential benefits of placing the canard in a non-coplanar (either high- or 
low-canard) configuration, a more detailed analysis of the lift and drag results is required. Polynomial 
curve fits through the computational and experimental results are used to determine the incremental 
lift and drag differences due to vertical positioning of the canard. 

Figure 3.10(a) shows the effect of canard vertical position on the incremental change in lift coefficient 
(A Cl) at different angles of attack. The reference line (A C L = 0.0) is the C L corresponding to either 
computations or experiment for the mid-canard configuration. Although the absolute A Cl values are 
not exactly predicted by the computations, it is noted that relatively small values are being compared. 
The qualitative trends indicating the relative effects of canard vertical position are captured well by the 
computations. Both the computations and experiment show increased lift at higher angles of attack 
for the high-canard case. However, the low-canard A Cl curves show significant loss of lift as a is 
increased. 

The incremental changes in drag coefficients (A C d ) given in figure 3.10(b) show that the high- 
canard configuration exhibits the best overall drag characteristics. At a higher Cl, the high-canard 
case has a lower Cd than both the mid- and low-canard cases. The results in figure 3.10 show clear 
evidence of an unfavorable canard-wing interaction for the low-canard configuration. 

3.2.2 Canard deflection angle 

Comparisons of the computed lift, drag and pitching moment curves with experimental data (ref. 12) 
are given in figure 3.11 for mid-canard deflections (<5 C ) of 0 and 10 deg. Figure 3.11(a) shows the 
nonlinearity in the lift curves. Note that the subtle differences in the experimental 6 C = 0 deg and 
8 C = 10 deg lift curves are captured well by the computations. The drag polar curves in figure 3.11(b) 
show that the undeflected mid-canard (<5 C = 0 deg) exhibits a lower overall drag coefficient. The drag 
polars also show that minimum Cd for the 10 deg case occurs at a negative lift coefficient. Figure 3.11(c) 
shows that increasing mid-canard deflection increases overall pitching moment for all computed angles 
of attack. Significant nonlinearities due to the canard-wing interaction are also noted in the pitching 
moment curves. 

Similar to the vertical-offset studies, the lift and pitching moment curves are separated into two 
configuration component regions in order to examine the effects of mid-canard deflections in more 
detail. Figures 3.12(a) and (b) illustrate lifts and moments, respectively, for the 6 C = 10 deg case. 
The wing portion lift curve of figure 3.12(a) shows the effect of the canard downwash on reducing 
wing lift. For example, at a = 0 deg, the additional canard lift (due to its deflection) is mostly 
cancelled by the negative lift generated by the wing (due to the canard downwash). The pitching 
moment curves in figure 3.12(b) show that the wing, in the presence of the canard, exhibits statically 
stable characteristics while the canard pitching moment causes the total configuration to be statically 
unstable. Such characteristics are often typical of canard-configured aircraft. Favorable comparisons 
in figures 3.12(a) and (b) between computations and experimental data indicate that the distribution 
of lift and moment between the component regions of the mid-canard deflection cases are captured well 
by the computations. 

An analysis of the effect of canard deflection on lift, pitching moment, and, particularly, trimmed 
conditions is of interest. For long-coupled canard configurations, the effect on aerodynamic performance 
is usually limited to the lift of the canard itself and the total pitching moment. However, the close- 
coupled canard has significant influences on wing performance as well. 

Figure 3.13 shows lift and moment curves for the canard and wing portions at various canard 
deflection angles ( 6 C ) ranging from -10 to +10 deg. While the canard portion lift curve (fig. 3.13(a)) 
exhibits increasing lift for increasing 8 C at a given a. the wing portion lift curve (fig. 3.13(b)) decreases 
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for increasing 6 C . As 8 C is increased, the canard generates greater lift and, consequently, greater 
downwash in the canard wake. The resulting effective wing angle of attack (a e ff) on the inboard 
portion of the wing decreases with increasing 6 C . Vortex visualization results presented later in this 
study will show that, for certain a and 6 C values, a e jj can be negative in the inboard wing while 
freestream a is positive. 

The canard and wing portion pitching moment curves for various canard deflection angles are given 
in figures 3.13(c) and (d), respectively. The canard portion pitching moment increases as 6 C is increased 
while the wing portion pitching moment is relatively insensitive to 8 C . The strong 6 C coupling with wing 
portion lift (fig. 3.13(b)) and the weak coupling with wing moment (fig. 3.13(d)) indicate a significant 
shift in wing center-of-pressure as a function of 6 C . Due to the canard downwash effect of “unloading'* 
the inboard portion of the wing, the aft-swept wing’s center-of-pressure location moves outboard and 
aft for increasing canard deflections. 

Figure 3.14 illustrates lift and moment as a function of 8 C at a representative a of 4.27 deg. The 
lift curves in figure 3.14(a) are given for the total configuration, canard portion and wing portion. For 
increasing 8 C , the cancellation of the lift between the canard and wing portions is clearly observed. 
At this moderate angle of attack, the total lift is minimally affected by the deflection of the canard. 
The moment curves in figure 3.14(b) show that the canard portion moment is the major contributor 
to changes in the total configuration moment. At this condition ( M,^ = 0.85, a = 4.27 deg), the 
configuration is trimmed about the moment center (69.17 cm from the nose) at 8 C —4 deg. 

For a realistic aircraft configuration, the determination of trimmed lift and drag is critical in assess- 
ing aerodynamic performance. A contour plot of pitching moment about the model center-of-gravitv 
location (c.g. = 59.14 cm from the nose) is given in figure 3.15. The trimmed (C m — 0.0) pitching mo- 
ment curve is highlighted and indicates that relatively small, negative canard deflections are required 
to trim this configuration at moderate angles of attack. 


3.3 Canard-Wing Vortex Interaction 

For the configurations in this study, many of the complex non-linear effects of the various canards on 
aerodynamic performance can be directly attributed to the canard and wing vortical flowfields and. in 
particular, the canard-wing vortex interaction. The current configuration flowfield is characterized by 
strong leading-edge vortices from the canard and the wing. As vortex strength increases, secondary 
separation induced by adverse pressure gradients can result in strong secondary vortices. The presence 
of secondary vortices will also serve to alter the position and trajectory of the primary canard and/or 
wing vortex. In cases of high vortex strength, tertiary separations can be observed as well. 

One of the primary mechanisms for the canard- wing interaction is the canard's influence on the wing 
leading-edge vortex. Both the canard downwash and the canard leading-edge vortex have pronounced 
effects on the formation arid subsequent trajectory of the wing vortex. The canard vortex influences 
the surface pressures on the wing directly and, to a greater extent, indirectly through its influences 
on the wing vortex. Furthermore, the effect of the canard vortex on the wing is highly sensitive to 
vortex strength and relative location, which are both significantly modified by variations in the canard 
position and deflection angles. 

3.3.1 Coplanar canard 

A baseline analysis of the canard-wing vortex interaction is performed for the undeflected mid-canard 
(coplanar) configuration at Moo = 0.90, Rec = 1.52 million, and various angles of attack. Comparisons 
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are made between canard-on and canard-off cases to determine the effects of the canard on the structure 
of the wing flowfield. 

Figure 3.16 illustrates a comparison between the upper-surface pressure contours of the canard- 
off and canard-on configurations at three angles of attack. These contours show a significant canard 
influence on the body as well as the influence of the canard and wing vortices on surface pressures. 
The vortex-induced low pressure region (L) gives a qualitative estimate of the vortex trajectory over 
the upper surface. As angle of attack increases, these low pressure regions move inboard. Note that, 
since the canard sweep angle is lower, and because there is an absence of a downwash field upstream 
of the canard, the recovery shock is strongest and most easily observed on the canard at the higher 
angles of attack. 

Although evident in the surface pressure contours of figure 3.16, the canard’s effect on the formation 
of the wing vortex is more clearly observed in the computed surface flow patterns given in figure 3.17. 
The secondary separation line (£> 2 ) due to the wing leading-edge vortex is observed for both the canard- 
on and canard-off cases. As the angle of attack is increased, the secondary separation and primary 
attachment lines move inboard and are indicative of the upward and inward movement of the wing 
primary vortex. The surface flow patterns for the canard-on cases show the modified wing flowfield 
including the delay in wing primary vortex formation. At a « 4 deg, no evidence of a leading-edge 
vortex is detected on the inboard portion of the wing. At a « 8 deg and 12 deg, the corresponding 
inboard surface flow patterns are influenced by a small leading-edge vortex which is comparable to 
a leading-edge vortex formed for the canard-off case at lower angles of attack. These observations 
correlate with the surface pressure distributions which were given earlier in figures 3.1 and 3.7. 

The dominating effect of the canard downwash on the wing, inboard of the canard-tip span-line, 
has already been illustrated. Outboard of the canard-tip span-line, the canard leading-edge vortex 
is the primary mechanism for the canard’s influence on the wing flowfield. Crossflow-plane vortex 
visualizations in the form of normalized total pressure contours and velocity vectors at two constant 
streamwise stations are given in figure 3.18, and show the effect of the canard-vortex induced flow on 
the wing. The approximate total pressure loss at each canard or wing primary vortex core is also given 
in figure 3.18 and is an indication of the instantaneous vortex strength. Although the computed canard 
vortex (fig. 3.18(a)) is considerably weaker than the corresponding wing vortex, its influence on the wing 
flowfield is still significant. For the two co- rotating vortices, each with a counter-clockwise rotation and 
position as shown in figure 3.18(a). the wing vortex flowfield induces a relative downward and inward 
motion of the canard vortex while the canard vortex induces an upward and outward movement of 
the wing vortex. In the absence of such interaction, the canard and wing vortex trajectories would be 
expected to follow an upward and outward path which would be dependent on the angle of attack and 
respective sweep angles. 

In addition to the canard-wing vortex interaction, the flowfield is further complicated by the pres- 
ence of a counter-rotating secondary vortex. The secondary vortex for the canard-off case is clearly 
visible in both total pressure contours and velocity vectors of figure 3.18(b). The smaller wing vortex 
for the canard-on case renders the corresponding secondary wing vortex undetectable in the scale of 
figure 3.18(a). However, computed surface flow patterns given earlier in figure 3.17 clearly show the 
existence of secondarv separation for both the canard-on as well as the canard-off cases. Evidence of a 
tertiary wing or canard vortex is not observed. 


3.3.2 Canard vertical position 

For the same flow conditions as the mid-canard cases given in the previous section (3.3.1), the differences 
in wing surface pressures between the high- and low-canard configurations are shown in figures 3.19(a) 
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and (b) for a = 4.27 deg and 8.55 deg. respectively. The upper and lower surface pressure distributions 
are given at three wing span stations for each angle of attack. All three stations are inboard of the 
canard-tip span line and, therefore, are directly influenced by the high- or low-canard downwash. 
Results in figure 3.19 show that the relative canard and wing positions have a large effect on the level 
of canard- wing interaction. 

At a = 4.27 deg (fig. 3.19(a)). small differences in wing surface pressures are observed. For the 
high-canard case, lower pressures on the upper surface near the wing leading edge indicate the possible 
formation of a leading-edge vortex. In the low-canard case, formation of the wing leading-edge vortex 
is inhibited due to the canard downwash and the relative positions of the canard and wing. 

In figure 3.19(b), the wing leading-edge vortex is clearly visible and contributes substantially to 
the lift of the high-canard configuration. However, even at this higher angle of attack, the low-canard 
results do not show evidence of a strong vortex on the inboard wing. In addition to inhibiting the 
formation of this vortex, the low-canard case also shows a substantial low pressure region on the wing 
lower surface. Both of these factors (absence of the wing vortex and low pressure on the wing lower 
surface) contribute to the loss of lift for the low-canard case observed in figures 3.8(a) and 3.10(a). 

Further insight into the unfavorable nature of the low-canard configuration can be gained by ex- 
amining the canard-wing flowfield structure. Figure 3.20 illustrates crossflow contours of stagnation 
pressures for the high- and low-canard configurations at a = 4.27 deg. For each configuration, five 
streamwise stations are chosen to illustrate the development and convection of the canard and wing 
vortices. The high-canard results (fig. 3.20(a)) show the formation of the canard vortex and its sub- 
sequent trajectory over the wing surface. At x/c = 3.2 and 3.6, formation of the wing leading-edge 
vortex is observed. Although the vortex positions differ, the high-canard-wing vortex interaction has 
similar characteristics to that of the mid-canard configuration shown in section 3.3.1. 

The low-canard case given in figure 3.20(b) shows a significantly different flowfield structure. Like 
the high-canard case, the low-canard vortex forms as visualized at x/c — 2.0. However, due to the 
relative position of the low-canard to the wing and the canard-wing aerodynamic interaction, the 
subsequent canard vortex is convected under the wing surface. The low-canard vortex can therefore 
induce an unfavorable low pressure region on the wing lower surface. The formation of the wing vortex 
on the upper surface is again observed in figure 3.20(b) at x/c = 3.6. 

Figure 3.21 shows that at the higher angle of attack of 8.55 deg. the low-canard vortex interacts 
directly with the wing. The original low-canard vortex (x/c = 2.0) splits into two (upper and lower 
canard) vortices upon impact with the wing leading edge (x/c — 2.8). Further downstream at x/c > 3.2. 
the split canard vortices continue to have the same rotational structure, but now exhibit independent 
trajectories. The canard lower vortex (ClL), with its counter-clockwise rotational sense as viewed in 
the figure, tends to move inboard due to its proximity to the wing lower surface. The canard upper 
vortex (C1U), also with a counter-clockwise rotation, is additionally influenced and moved inward by 
the formation of the wing vortex (Wl) on the wing upper surface (x/c = 3.2 and 3.6). 

To visualize the results in figure 3.21 more clearly, a qualitative schematic of the canard- wing vortex 
trajectories for the low-canard case is given in figure 3.22. The low-canard vortex impacts the wing 
leading-edge and is split into the upper and lower vortices. The wing vortex forms on the wing upper 
surface and interacts with the canard upper vortex. 

A perspective view of the low-canard vortex interaction with the wing is given in figure 3.23 for 
a — 8.55 deg. Crossflow contours of stagnation pressures show the formation of the canard vortex, 
the splitting of the canard vortex by the wing, and the formation of the wing vortex on the outboard 
wing. The effect of this wing vortex on the canard upper-vortex trajectory is also clearly visible. 
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3.3.3 Canard deflection angle 

In addition to those canard influences due to changing vertical position, the additional freedom-of- 
motion of mid-canard deflections allows for further variations in the types of canard-wing vortex struc- 
tures. Changing 6 C affects the strength of the canard downwash and vortex, the location of the canard 
trailing-edge relative to the wing, and the subsequent canard vortex trajectory in the wing flowfield. 

Results given in section 3.2.2 illustrated the effect of canard deflections on lift, drag and pitching 
moments for the entire configuration and the separate canard and wing component regions. The 
aerodynamic performance characteristics from that section are better understood with further analyses 
of the canard-wing vortex interactions. The differences in surface flow patterns for 8 C = 0 deg, 

5 deg and 10 deg are illustrated in figure 3.24. As for the undeflected mid-canard case, the secondary 
separation lines from the canard and wing primary vortices are indicated (S2) for each case and show 
that viscous effects are significant to the overall flow structure. As canard deflection is increased, the 
secondary separation line on the canard moves inward reflecting the upward and inward shift of the 
canard vortex. In all three cases, the secondary separation line forms near the apex of the canard while 
its formation on the wing inboard is inhibited. At this angle of attack, the formation of the wing vortex 
is delayed until outboard of the canard-tip span location. 

Figure 3.25 illustrates the effect of canard deflection on wing surface pressure at 35, 55 and 75 per- 
cent semispan. For clarity, the upper surface pressure curves are indicated (U). As canard deflection is 
increased from 0 to 10 deg, the near-leading-edge lower surface pressure at 35 percent span transitions 
from positive to negative values. At — 10 deg, a strong suction peak exists on the lower surface, 
indicating the possibility of a wing vortex on the wing lower surface. Furthermore, high upper surface 
pressures at 35 percent span are clearly indicative of a decrease in local angle of attack for the inboard 
wing with increasing canard deflection angle. The loss of wing lift due to these pressure changes was 
shown previously' in figure 3.13(b). Towards the outboard stations, figure 3.25 shows that the effect of 
canard deflection on surface pressure is minimal. 

To visualize the effect of canard deflection on the canard-wing vortex interaction, crossflow contours 
of total stagnation pressures are given in figure 3.26 for 8 C = 5 deg and 10 deg. 1 he crossflow contours 
are shown at x/c = 2.0. 2.8 and 3.5 as measured from the fuselage nose. The relative primary vortex 
strengths are indicated in the figure by the approximate total pressure loss in the vortex core. For 
S c = 5 deg (fig. 3.26(a)), there are canard primary (Cl), canard secondary (C2), and wing primary 
(W1U) vortices over the upper surfaces of the canard and wing. The distortion of the canard primary 
vortex at x/c = 2.8 is due. in part, to the canard secondary vortex. 

At 8 C = 10 deg (fig. 3.26(b)). an additional wing primary vortex (W1L) forms and is visible on 
tlie wing lower surface (x/c — 2.8). This vortex causes the suction peak on the lower surface of the 
wing observed earlier in figure 3.25 and would generally be considered an unfavorable effect. At these 
flow conditions, there is a transition from negative to positive effective wing angle of attack, a e //, from 
the wing-root to wing-tip locations. At the last station, x/c = 3.5, the wing flowfield structure for the 
8 C = 10 deg case includes two wing primary vortices (upper and lower) as well as the canard primary 
vortex. 

The mutual interaction between these vortices is an important factor in wing performance and is 
complex in nature. A comparison of the relative position of the canard vortex between x/c = 2.8 and 
3.5 in figure 3.26 shows that the trajectory of the canard vortex is influenced by the wing upper-surface 
vortex. For two co- rotating vortices (canard and wing primary), each with a position and counter- 
clockwise rotation as shown in the figure, the wing vortex induces a relative inward motion of the 
canard vortex. 
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The influence of the wing vortex on the canard vortex trajectory is more clearly visualized in 
figure 3.27, which illustrates the upper surface canard and wing vortex systems for 6 C — 5 deg and 
10 deg. Crossflow contours of total stagnation pressure are shown at selected streamwise stations. As 
shown earlier in figures 3.24 to 3.26 for both 6 C = 5 deg and 10 deg, the formation of the wing vortex 
on the wing upper surface is inhibited by the canard downwash until outboard of the canard-tip span 
location. The difference in the canard primary and secondary vortex strengths between the two 6 C cases 
is qualitatively observed. In each case, as the canard primary vortex is convected over the w'ing upper 
surface, the wing primary vortex in the outboard region clearly modifies the canard vortex trajectory 
inward. 

Figure 3.28 shows a perspective view of the overall flowfield for the configuration at a = 4.27 deg 
and 6 C = 10 deg. The darker region over the canard upper surface indicates low pressure induced by 
the canard vortex. The particle traces released near the wing root and towards the wing tip show the 
formation of the leading-edge vortices on both the upper and lower surfaces. 

3.4 Effect of Canard on Wing Vortex Breakdown 

The potential for a canard to delay or eliminate wing vortex breakdown is of significant interest and has 
been the topic of numerous experimental studies (refs. 19, 21, 53, and 54) and a recent numerical study 
(ref. 55). By comparing vortex lift theory with experimental results, Gloss and Washburn (ref. 19) 
found that wing vortex burst occurs at a % 13 deg for the current canard-off case. For the undeflected 
mid-canard case, their study indicated no evidence of wing vortex burst for angles of attack up to at 
least 20 deg. 

Figure 3.29 shows a crossflow-plane visualization of the predicted wing vortex for the canard-off 
case at three angles of attack. At the higher angles of attack, the resulting increased vortex strength 
and relative shift of the core location is observed in the normalized total pressure contours. The 
development of the secondary vortex is also noted. The corresponding off-surface particle traces for the 
canard-off cases are shown in figure 3.30. Evidence of vortex burst in the computational simulation is 
observed at a = 12.38 deg. Crossflow- plane visualization of scaled axial velocity contours in figure 3.31 
serves to confirm the presence of predicted vortex breakdown over the wing for the canard-off case. 
Reversed axial flow in the core of the primary vortex is identified (ref. 56) as a qualitative indication 
of vortex breakdown, and is observed in figure 3.31. From figure 3.30, the wing vortex burst location 
appears to be near the trailing edge of the wing and indicates that the computed angle of attack for 
vortex burst is approximately 12 deg, which is within reasonable agreement with the experimental 
observations (ref. 19). 

Computed particle traces for the undeflected mid-canard case at a = 12.38 deg are given in 

figure 3.32. The lower leading-edge sweep-angle of the canard results in a canard vortex burst further 
upstream of the canard trailing edge compared to the wing vortex burst relative to its trailing edge 
shown earlier. Figure 3.32 shows a stable wing vortex with no evidence of wing vortex breakdown in 
the presence of the undeflected mid-canard. Although not computed directly for this study, the results 
shown in Secs. 3.2 and 3.3 indicate that wing vortex breakdown is quite possible given certain canard 
vertical positions and deflection angles. 

The good comparisons with the experimental data for both the canard-on and canard-off cases at 
q ^ 12 deg and the corresponding canard-off vortex burst location near the wing trailing edge show 
that the computed relative effects on vortex breakdown location and behavior in the presence of the 
canard are accurate. However, it is noted that vortex breakdown is generally an unsteady phenomenon 
and steady-state computations downstream of the burst location may be suspect. To perform a detailed 
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Figure 3.1. Comparison of computed and experimental steady-state surface pressure co 
efficients for the baseline grid with and without mid-canard. A/oo = 0.90, a ~ 4 

Re £ = 1-52 million. 
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Figure 3.2. Comparison of computed and experimental force coefficients for the baseline 
zrid with and without mid-canard. ~ 0.90, Re = = 1.52 million. 
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Figure 3.3. Comparison of component force coefficients for the baseline grid with and 
without mid-canard. (Lift and moment curves are given for shaded regions of the geometry). 
Mo© - 0.90, # e c = 1-52 million. 
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Figure 3.4. Comparison of baseline and refined grid surface pressure coefficients with ex- 
periment. M a o = 0.90, a = 8.21°, Re d = 1.52 million. 
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Figure 3.5. Comparison of baseline and refined grid surface pressure coefficients with ex- 
periment. Moo = 0.90, a - 12.38°, Re d = 1.52 million. 
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Figure 3.6. Comparison of baseline and refined grid force coefficients with experiment 
M~, = 0.90. Rer = 1.52 million. 
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Figure 3.8. Comparison of computational and experimental forces and moments for the 
high-, mid- and low-canard configurations. = 0.90. Rec = 1*52 million. 
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Figure 3.9. Comparison, of component lift and pitching moment curves for the high-canard 
configuration. M ^ — 0.90, Re f = 1.52 million. 
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Figure 3.10. Effect of canard vertical position on the incremental change in lift and drag 
coefficients. = 0.90, Re £ = 1.52 million. 
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Figure 3,11. Comparison of computed and experimental forces and moments for the un- 
deflected ( 6 C — 0,0°) and deflected (<5 C = 10.0°) canard configurations. = 0.85, 

Re d ~ 2.82 million. 
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Figure 3.12. Comparison of component lift and pitching moment curves for the deflected 
canard configuration at 6 C = 10.0°. Mcc = 0.85, Re £ = 2.82 million. 


51 





Canard Deflection, 5 C 
Computations Experiment 


10 degrees □ 




a) Canard portion lift c) Canard portion moment 




Figure 3.13. Effect of canard deflection, <5 C , on wing and canard portion lift and moment 
curves. = 0.85. Re e = 2.82 million. 
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Figure 3.14. Effect of canard deflection, <5 C , on component lifts and pitching moments. 
M 0 o = 0.85, a = 4.27°, Re s = 2.82 million. 
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Figure 3.15. Contour plot of computed total pitching moment with moment center at c.g. 
location (59.14 cm from nose). M 0 0 = 0.85, Re 5 = 2.82 million. 
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Figure 3.16. Computed upper surface pressure contours for the wing-body configuration 
with and without mid-canard. = 0.90, Re s = 1.52 million. 
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a) Canard-on b) Canard-off 

Figure 3.18. Crossflow plane visualization of normalized total pressure contours and scaled 
velocity vectors of the canard and wing leading-edge vortices. - 0.90, a = 8°, 

Re c = 1.52 million. 
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Figure 3.20. Crossflow plane visualization of normalized total pressure contours for high- 
and low-canard configurations. = 0.90, a — 4.27°, Re 5 = 1.52 million. 
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Figure 3.21. Crossflow plane visualization of normalized total pressure contours for the 
low-canard configuration. Moo = 0.90, a = 8.55°, Re ? - 1-52 million. 
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Figure 3.22. Schematic of the canard and wing vortex trajectories for the low-canard con- 
figuration. 



Figure 3.23. Perspective view of normalized total pressure contours for the low-canard 
configuration. = 0.90, a = 8.55°, Re 5 = 1.52 million. 
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Figure 3.26. Crossflow plane visualization of normalized total pressure contours for canard 
deflections of <5 C = 5° and 10°. Moo = 0.85. a = 4.27 Re- C = 2.82 million. 
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Figure 3.27. Perspective views of upper surface canard-wing vortex structure for canard 
deflections of 6 C - 5° and 10°. M ^ = 0.85. a - 4.27°, Re? = 2.82 million. 
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Figure 3.29. Crossflow plane visualization of normalized total pressure contours of the wing 
leading-edge vortex for the canard-off case. A/^ = 0.90. Re ? = 1.52 million, x/c — 3.6. 
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Figure 3.31. Crossflow plane visualization of scaled axial velocity contours for the canard-off 
case. Moo = 0.90. Re g = 1-52 million, x/c = 3.6. 




Figure 3.32. Off-surface particle traces of the canard and wing vortices for the canard-on 
case. Moo = 0.90, a = 12.38°, Re? = 1.52 million. 
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Chapter 4 

Effect of Fixed Canard on Unsteady 
Aerodynamics 


The effects of fixed canards on the unsteady longitudinal aerodynamics of the canard-wing-body con- 
figuration are presented in this section. As outlined earlier in the Introduction, many benefits of 
canards are realized during maneuver and other dynamic conditions which result in highl\ unsteady 
aerodynamics. In such conditions, the fixed canard has a pronounced influence on the unsteady forces, 
moments, and ultimately the performance and stability of the configuration. 

Computational simulations solving the time-accurate Reynolds-averaged Navier-Stokes equations 
are used to study the effects of the canard on the pitch-up ramp motion and pitch oscillations of 
the canard- wing-bodv configuration. All fixed canard simulations in this section are performed for 
the undeflected mid-canard configuration. Comparisons are made with available experimental data 
(ref. 24). and numerical convergence is demonstrated through spatial and temporal refinement studies. 

The results of the simulations are used to study the unsteady forces and moments due to the canard, 
and the unsteady canard-wing vortical flowfield. Ultimately, in section 4.3, the capability of utilizing 
high-order CFD for configuration stability evaluation is demonstrated through a dynamic stability 
analysis (derived in sec. 2.7) of the time-accurate pitch-oscillation predictions. 

4.1 Pitch-Up Ramp Motion 

Results are presented to provide the analyses of the unsteady canard-wing-body flowfield associated 
with pitch-up ramp motions. All unsteady results in this section are computed at M 0 0 = 0.90 and a 
Reynolds number based on the mean aerodynamic chord of the wing (Reg) of 1.52 million. Convergence 
of the unsteady flowfield is verified using both spatial and time-step refinement studies. 

All unsteady canard- wing-body results are computed for a from 0 deg to either 12.83 deg or 15 deg 
and for non-dimensional pitch rates (defined by A p — ac/Li x ) of 0.05 and 0.10. Ramp motions are 
started from converged steady-state solutions at an initial angle of attack (a*) and held at the final 
angle (ay) for a fixed length of time. The pitch axis of the ramp motion and the pitching moment 
results are taken from the model c.g. location shown in figure 1.4. 

A typical ramp motion from a* = 0 deg to ay = 15 deg is illustrated in figure 4.1. In order 
to directly compare unsteady results at different pitch rates and steady-state results, time (t) is given 
in degrees. During the ramp motion (0 deg < t < 15 deg) a and t are equal (a = t). However, for 
t > 15 deg, a is held constant at 15 deg (a = 15 deg or ay). It is noted that during the ramp motion. 
a increases linearly as a function of time. Consequently, the ramp motion has an impulsive start and 
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stop which significantly influences the unsteady results. (Note that the conversion to physical time can 
be easily obtained using the definition A p = ac/Ucc and replacing a with At ra d/ At sec .) 

4.1.1 Aerodynamic performance 

Computed time histories of lift, drag and pitching moment for the configuration, with and without 
canard, undergoing a ramp motion (pitch rate, A p = 0.10) are illustrated in figure 4.2. Computational 
and experimental steady-state results are also given for the corresponding angles of attack. The lift 
curves of figure 4.2 show a significant lift increase for the unsteady canard-on case. Early in the motion 
(t < 4 deg), the unsteady canard-on lift is slightly lower than the canard-off lift due to the pitch-rate 
effects of the configuration. Since the entire canard is forward of the pitch axis, there is an initial loss 
of lift on the canard at the start of the ramp motion (t = 0 deg). As the ramp motion continues, the 
canard-on case exhibits increased lift over both the canard-off and steady-state canard-on cases. 

The time histories of drag coefficient show that at a given instantaneous angle of attack, there 
is higher drag for the canard-on case. However, by replotting the drag results of figure 4.2 in drag- 
polar form, figure 4.3 shows that the unsteady canard-on case exhibits improved dynamic lift-to-drag 
performance at the higher angles of attack. 

The pitching moments given in figure 4.2 also illustrate the significant influence of the fluid virtual 
(or apparent (ref. 57)) mass at the .beginning (t = 0 deg) and end (t = 15 deg) of the ramp motion. 
The virtual mass acts to counter the acceleration of the body and, therefore, causes a rapid nose-down 
pitching moment at t = 0 deg and a nose-up moment at t = 15 deg. Beyond t = 15 deg, the lift, 
drag and pitching moment values converge toward the steady-state a — 15 deg result. 

The effect of pitch rate on lift for the canard-on case is given in figure 4.4 for A p = 0.10 and 0.05. 
Note that since comparisons are made at instantaneous angles of attack during the ramp motion (t, in 
deg), the physical time (t se c) between the two pitch rates differ by a factor of two. Figure 4.4 shows 
that the lift of the canard configuration is increased at the higher pitch rate throughout the ramp 
motion. Further into the ramp motion ( t > 10 deg), the pitch rate and total lift are approximately 
proportional at given angles of attack. 

Better insight into the dynamic loads produced by the ramp motion can be attained by examining 
the separate component regions of the geometry. Similar to the steady-state component results given 
earlier, the canard region consists of the canard and the body forward of the wing leading-edge root 
location (fore-body). The wing region consists of the wing and the remaining aft-body (not including 
the sting). 

Time histories of the canard and wing region lift contributions for the A p = 0.05 ramp motion 
are illustrated in figure 4.5. Computed steady-state lift coefficients are also given for reference. By 
definition, the total configuration lift is the sum of the canard and wing region lift. The effects of the 
fluid virtual mass are evident at both a t = 0 deg and aj = 15 deg. Due to the relative locations of 
the canard and wing to the pitch axis, there is an initial increase in wing region lift and a decrease in 
canard region lift at cq = 0 deg. These trends are then reversed at ocf = 15 deg. 

Figure 4.5 shows that, in comparison to the steady-state results, there is a net loss of lift for the 
canard throughout the ramp motion. The increased lift for the total configuration during the ramp 
motion is due to the large increase in the lift of the wing. Although the canard-wing flowfield is highly 
non-linear and complex in nature, it can be quite useful to qualitatively describe the effects of the 
unsteady motion in terms of simplified linear theory. In linear theory, the lift of the wing or canard can 
be written as the sum of the static lift (Cx a a), the lift due to pitch or quasi-static lift ( Ci q q ), and the 
unsteady lift (C^ q q). The terms Cl q , Cx 9 , and Cl 6 are the stability derivatives (dCi/da) o. (dCi/dq ) o 
and (dCi/da) o, respectively; q is the angular velocity of the configuration about the pitch axis. The 
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notation () 0 indicates that partial derivatives are evaluated assuming the disturbance quantities are 
zero. Note that for the ramp motion, q = a. 

For the canard or wing, the lift contribution due to pitch is represented in the C Lq and C Ln terms. 
The loss of canard lift is primarily due to the canard location forward of the pitch axis which, if 
described by linear stability theory, causes a large negative contribution to CL canard from the C Lq g 
term for high pitch rates. The directions of the induced velocities due to pitch are given schematically 
in figure 4.5 as well. After the ramp motion stops ( t > 15 deg), the wing region lift contribution 
converges to the steady-state result much more slowly than the canard region lift. Results shown in 
section 4.1.3 indicate that the canard-wing vortex interaction is a significant factor in this phenomenon. 


4.1.2 Spatial and temporal convergence 

To confirm the convergence and accuracy of the pitch-up ramp motion computations, time-step and 
grid refinements are performed for the A v = 0.05 ramp case presented earlier in figure 4.5. Figures 4.6 
and 4.7 show the effect of time-step size (At in deg) on the unsteady lift curves near t = 0 deg and 
t = 15 deg. respectively, for the component regions of the geometry. In the previous computations 
presented in this section, At = 0.0025 deg was used. 

From figure 4.6, it is clear that with decreasing time-step size, the time-accurate solutions converge 
quickly. In fact, even when using a larger step size, the lift curves compare favorably for t > 0.1 deg. 
Figure 4.7 shows that, after the ramp motion has stopped, the transient solution is adequately converged 
in time using At = 0.0025 deg. 

Spatial accuracv is demonstrated in figure 4.8 by the favorable comparison between baseline and 
refined grid solutions of component lift curves. The refined grid contains over 1.7 million points in 
the flowfield and was used extensively in Chapter 3 for the steady-state study. Figure 4.8 shows that, 
with the refined grid, slightly higher lift coefficients are predicted due to the improved resolution of the 
canard and wing leading-edge vortices. The remaining analyses of the unsteady canard-wing vortex 
interaction are performed using the refined grid. 

4.1.3 Canard- wing vortex interaction 

An analysis of the canard- wing vortex interaction provides further understanding of the complex flow- 
field that produces the unsteady load characteristics of the ramp motions shown earlier. Computations 
are performed for the canard-wing-body configuration undergoing a ramp motion of a, = 0.00 deg to 
aj = 12.83 deg and at a pitch rate of A p = 0.05. Comparisons are also made with the corresponding 
steady-state results. 

The steady and unsteady ramp motion results are summarized in figures 4.9 and 4.10, respectively, 
which show instantaneous streamlines released on the upper surfaces of the canard and wing. Positions 
A s . B s . and C s , as indicated in figure 4.9, correspond to steady-state a = 4.21 deg. 8.55 deg, and 
12.83 deg, respectively. Similarly, during the ramp motion (fig. 4.10), positions .4, B, and C correspond 
to instantaneous a = 4.21 deg. 8.55 deg, and 12.83 deg. The ramp motion stops at a f = 12.83 deg 
and the configuration is held fixed for all remaining time (t > 12.83 deg). Time-accurate transient 
results are given for t > 12.83 deg at instances D, E, and F. The a vs. time plot given in figure 4.10 
illustrates the relative time instances of the ramp motion ( A-F ). Corresponding upper surface pressure 
contours for the steady and unsteady cases are given in figure 4.11 for the same instances during the 
ramp motion as in figures 4.9 and 4.10 (j4 s — C s and A — C, respectively). 

A comparison between the steady and unsteady results shows significant differences in the char- 
acteristics of the canard and wing vortices during the early portion of the ramp motion. Steady and 
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unsteady case total pressure contours at a crossflow plane of x/c = 3.6, as measured from the 

body nose, are given in figure 4.12 for these same angles of attack. As observed in figures 4.10, 4.11, 
and 4.12(b), a delay in canard vortex formation occurs for the unsteady case which correlates with the 
loss of canard lift due to the pitching motion (fig. 4.5). For the steady case, figure 4.12(a) shows the 
presence of the convected canard vortex, even at the lower angles of attack. Differences between the 
steady and unsteady wing vortices are also noted in figure 4.12. During the ramp motion, the increased 
strength of the wing vortex (fig. 4.12(b)) contributes to the increased lift of the wing portion observed 
earlier in figure 4.5. 

For the current ramp motion, as the configuration reaches a/, the motion stops. Figure 4.10 shows 
the formation of the canard vortex and, at later times, evidence of canard vortex breakdown. Total 
pressure contours in figure 4.13 illustrate the transient vortical flowfield of the canard (x/c = 2.0) and 
wing (x/c = 3.6) once the ramp motion has stopped. The development of the canard vortex, which 
was inhibited during the early portion of the ramp motion, is clearly observed. The canard vortex 
convects downstream and the changes in canard lift and downwash influence the characteristics of the 
wing vortex. As time increases from instances D to F (fig. 4.13), the canard- wing flowfield approaches 
the steady-state results illustrated earlier in figure 4.12(a). 

Although the canard vortex forms rapidly once the ramp motion has stopped, figure 4.13 shows that 
the increasing strength of the canard vortex, its subsequent convection downstream, and its interaction 
with the wing vortex occur more slowly. As also observed earlier in figure 4.5, it is evident that 
the time scale is considerably larger for the wing than the canard to reach a relatively steady-state 
condition. Although some unsteadiness remains in the flow due to the breakdown of the canard vortex, 
the time-accurate results indicate a relatively stable flowfield for large time (t >> 12.83 deg). 

4.2 Pitch Oscillation 

Results of pitch-oscillation simulations are presented to understand the unsteady flowfield, provide 
analyses of the canard-wing-body dynamic stability characteristics and to evaluate CFD accuracy and 
efficiency issues in predicting stability parameters. In order to match experimental test conditions 
(ref. 24) for the dynamic stability analysis, unsteady results are computed at Mqq = 0.70 with the 
Reynolds number based on the mean aerodynamic chord of the wing (Rec) of 3.05 million. 

Unsteady computations are performed for the mid-canard configuration undergoing pitch oscilla- 
tions about the pitch axis (shown earlier in fig. 1.4). Computations are made at mean angles of attack. 
a m . of 0 to 16 deg and a pitch amplitude, A, of 1.0 deg. A reduced frequency, as defined bv k = 7 /ff 

■ZC'oo 

of 0.236 radians is used in the majority of the computations. 

The computed time histories for the lift and moment coefficients are given in figures 4.14 and 4.15, 
respectively, as a function of non-dimensional time. ut. Each case was computed for three full cycles 
of oscillation about a mean angle of attack of 4 deg, w r ith the time histories of the latter two cycles 
shown in the figures. The angle of attack, with the appropriate scale on the right vertical axis, is given 
in each figure for reference. 

A comparison of the lift time histories in figure 4.14(a) show's that the canard-on and canard-off 
cases exhibit similar trends. The amplitude of the oscillation in lift is slightly higher and the average 
lift coefficient is slightly lowor for the canard-on case. Figure 4.14(a) also show's no significant phase 
shift for the total lift due to the canard. A breakdown! of the lift time histories into component regions 
is given in figure 4.14(b) for the canard-on case. As with the pitch- up ramp cases given earlier, the 
canard and wing regions include the fore- and aft-body portions, respectively. Figure 4.14(b) shows 
that significant phase differences exist between the canard and wing lift curves. In fact, the wing 
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lift is shown to lead the pitch motion while the canard lift appears to lag. Previous results given m 
section 4.1 demonstrated that the lead and lag in wing and canard lift, curves is primarily caused by 
the induced velocities due to the rapid pitch motion. A comparison of figures 4.14(a) (canard-off) and 
4.14(b) (wing lift) also shows that the presence of the canard causes a phase shift in the wing region 
lift curve. The resulting total lift for the canard-on case slightly leads the pitch oscillation. 

The time histories of canard-on/off total and canard-on component pitching moments are ghen in 
figures 4.15(a) and (b). respectively. Figure 4.15(a) shows that the canard-on case exhibits a higher 
average nose-up moment. In contrast to the total lift curves of figure 4.14(a). a more significant 
phase” shift between canard-on and canard-off total moment curves is also noted. The breakdown into 
components of the moment contributions (fig. 4.15(b)) illustrates the large phase differences between 
the canard and wing region moments. Due to the forward position of the canard, the nose-up pitching 
moment of the canard results in a significant nose-up shift of the total moment curve. 

In order to isolate the effects of the canard on the amplitude and phase of the lift/moment curves, 
a plot of the change in lift/moment as a function of the change in angle of attack is given in figure 4.16. 
In both the lift (fig. 4.16(a)) and moment (fig. 4.16(b)) curves, the change in lift/moment is taken from 
the average lift/moment coefficients for the respective canard-on or canard-off cases. The direction of 
the curves are also indicated in the figure. From figure 4.16(a), the larger oscillation amplitude and the 
minimal phase change of lift coefficient for the canard-on case is more clearly observed. Figure 4.16(b) 
shows the significant phase shift in moment coefficient for the canard-on case. 

Further insight into the energy balance related to the motion of the configuration can be gained 
by integrating the pitching moment curves of figure 4.16. For the co-located moment center and pitch 
axis, an evaluation of 

W p = J> C m (a)da ( 4 - 2 - 1) 

for a single cycle represents the amount of nondimensional work done by the fluid due to the motion 
of the configuration. A similar analysis for roll oscillations of delta wings has been performed in a 
previous computational study (ref. 58). For the pitching moment curves of figure 4.16(b), which have a 
counter-clockwise sense, an evaluation of equation 4.2-1 yields W p < 0 indicating that work (or energy) 
is imparted to the fluid by the pitch oscillation of the configuration. Therefore, the configuration is 
expected to be positively damped in pitch. Furthermore, the apparent increased area contained within 
the canard-on moment curve (fig. 4.16(b)) indicates that more work is imparted to the fluid for the 
canard-on case, leading to a higher damped condition. 

The periodicity of the oscillatory motion is demonstrated with the plot of lift and moment as a 
function of instantaneous angle of attack given in figure 4.17. The time-accurate computations are 
started from converged steady-state solutions at a m . Figure 4.17 shows that lift and moment time 
histories reach a periodic oscillatory response within 2-3 cycles. A Fourier analysis was also performed 
on the second and third cycles of the oscillatory motion and the results showed less than 0.5 percent 
difference between the cycles in computed Fourier coefficients. Therefore, for the current configuration 
and flow conditions, the second cycle of oscillation is sufficiently accurate to perform the detailed 
dynamic stability analysis. 

4.3 Dynamic Stability 

A comparison between the lift /moment curves and the appropriate Fourier series using the first three 
Fourier coefficients (commonly referred to as ao-fli- and bi) is given in figure 4.18 for the canard-on 
configuration. The coefficients were determined from a Fourier analysis of the computed second cycle 
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and the excellent agreement confirms that the higher-order terms can be neglected. These coefficients 
are then scaled with A and Ak as described in section 2.7 to obtain the corresponding dynamic stability 
parameters. 

Figure 4.19 shows a generally favorable comparison between the computed and experimental (ref. 24) 
dynamic stability parameters for the canard-off case at various mean angles of attack (a m ) and an 
amplitude, A = 1.0 deg. The stability parameters based on normal force are given in figures 4.19(a) 
and (b) as the pitch rate and pitch displacement parameters, respectively. Damping-in-pitch and 
pitch-oscillatory stability, which are based on pitching moment coefficients, are given in figures 4.19(c) 
and (d), respectively. In addition, results from fine grid computations at two different a m are given in 
figure 4.19 and show improvement in the damping-in-pitch parameter (fig. 4.19(c)) at a m ~ 16 deg. For 
moderate angles of attack (a < 10 deg), both computations and experiment predict positive damping, 
(Cm Q + C m ) < 0 in figure 4.19(c), and positive oscillatory stability, ( C ma ~~ k 2 C mq ) < 0 in figure 4.19(d), 
for the canard-off configuration. 

4.3.1 Effect of canard 

The canard is expected to have a pronounced effect on the dynamic stability of the wing-body con- 
figuration. Results from Chapters 3 and 4 on this configuration have shown that the close-coupled 
canard has a large influence on both the steady-state and unsteady wing aerodynamics. The effect 
of the mid-canard on longitudinal dynamic stability parameters is summarized in figure 4.20. Each 
data point on the curves represents the predictions at a mean angle of attack of a m and an amplitude, 
A = 1.0 deg. 

Figure 4.20 shows that the presence of the canard significantly changes the dynamic stability charac- 
teristics of the configuration. In particular, the pitch-rate parameter in figure 4.20(a) shows an increase 
throughout the mean angle-of-attack range, while the pitch-displacement parameter (fig. 4.20(b)) shows 
a dependence on angle of attack for the effect of the canard. For the pitching moment based parameters 
(figs. 4.20(c) and (d)), the computations show an increase in damping and a decrease in longitudinal 
stability with the canard. 

The changes in damping-in-pitch and oscillatory stability due to the presence of the mid-canard 
is of particular interest and is again summarized in figure 4.21. Differences between canard-on and 
canard-off stability parameters are shown in figure 4.21 with the “zero-line'' representing the canard-off 
case. As with figure 4.20 earlier, figure 4.21 also illustrates that the canard increases the damping-in- 
pitch of the configuration while reducing the oscillatory stability. The increase in damping with the 
canard correlates with the earlier analysis of the fluid “work’' (eq. 4.2-1) described in section 4.2 and 
shown in figure 4.16(b). 

For reference, the experimental data for the high-canard configuration is shown in figure 4.21. 
Since there is a difference in the canard vertical positions between the computation (mid-canard) 
and experiment (high-canard), the comparison with the experimental measurements is intended as a 
qualitative confirmation of the computational findings. However, it is noted that the agreement with 
experiment is quite favorable, except at a m = 0 deg, where the wake of the mid-canard in the 

computational model directly interacts with the co-planar wing. 

The results given above are based on the computations of the canard-off and canard-on config- 
urations undergoing pitch oscillations at a reduced frequency, k = 0.268. On the other hand, the 
experiment (ref. 24) was conducted at a one-order-of-magnitude lower reduced frequency, k — 0.0268. 
The implications of this difference and, in general, the effect of reduced frequency on predicting dynamic 
stability parameters is a topic of considerable interest and is discussed in the following section. 
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4.3.2 Reduced frequency 

Time-accurate Navier-Stokes computations are generally restricted to small time-step sizes by accuracy 
and numerical stability constraints. As a result, simulations of low reduced-frequency cases, which are 
typical of small-scale wind-tunnel tests, can require extensive amounts of computer resources. The 
number of iterations required for each cycle of oscillation becomes inversely proportional to reduced 
frequency. Therefore, in order to improve computational efficiency, it is highly desirable to evaluate 
the effects of reduced frequency on the accuracy of the dynamic stability parameter predictions. 

In particular, such an evaluation can provide guidance on the highest reduced frequencies which 
can be used in computational simulations for a given configuration and flow condition. When the 
simulated reduced frequencies become too high, the higher-order terms of equations 2.7-1 and 2.7-7 
become significant and can no longer be neglected. Furthermore, the classical assumption of linear 
aerodynamic theory from which equation 2.^-1 is based, can no longer be taken. 

Figure 4.22 illustrates the effect of reduced frequency (k) on the normal force and pitching moment 
coefficients as a function of angle of attack for the canard-on configuration. The mean angle of attack 
is 4 deg and reduced frequency varies from k = 0.268 to 1.072. The increase in amplitude and shift in 
phase due to increasing k is observed in figures 4.22(a) and (b) for both the normal force and pitching 
moment coefficients, respectively'. 

From the dynamic stability analysis, the effect of k on the damping-in-pitch and pitch oscillatory 
stability parameters is given in figure 4.23 for the canard-on case. For clarity, the reduced frequency 
(horizontal axis) is plotted as the inverse of k (1/k) on a log scale. The experimental values are given 
at the actual reduced frequencies (k) reported in the wind-tunnel tests. Figure 4.23 shows that the 
dynamic stability parameters are accurately computed for a wide range of reduced frequencies. In fact, 
for damping-in-pitch. a computed reduced frequency' which is up to 1.5 orders-of-magnitude higher 
than the experimental value continues to yield accurate stability' parameter predictions. In the present 
studv. a reduced frequency' of k 0.268 is utilized and allows for an order-of-magnitude reduction 
in the required computational time. For dynamic stability evaluations in general, a similar study on 
the sensitivities to reduced frequency has the potential to reduce computational expenses by orders of 
magnitude. 
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Figure 4.1. Schematic of unsteady pitch-up ramp motion from a z = 0° to aj = 15°. 
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Figure 4.2. Effect of canard on unsteady aerodynamic loads during pitch-up ramp motion. 
jV/oo = 0.90, A p = 0.10, cti = 0°, Q f = 15°, Re 5 = 1.52 million. 
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Figure 4.3. Drag polar of the canard-wing-body configuration undergoing pitch-up ramp 
motion. = 0.90, A p — 0.10, — 0°, a/ = 15°, Re d = 1.52 million. 
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Figure 4.4. Effect of pitch rate on canard- wing-body lift. = 0.90, a t = 0°, aj = 15°, 

Re d = 1.52 million. 
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Figure 4.5. Comparison of steady and unsteady component lift curves for the canard con 
figuration. M c 0 = 0.90. A p = 0.05, a, = 0°, a f = 15°, Re s = 1.52 million. 


At (deg); also Aa 
— 0.01000 

— 0.00500 

— 0.00250 



0 0.1 0.2 0.3 0.4 0.5 


time (deg) 


Figure 4.6. Effect of time-step size on the unsteady component lift curves near the beginning 
of the ramp motion, a t . Mqq - 0.90, A p = 0.05, a{ — 0°, aj = 15°, Re? = 1-52 million. 
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Figure 4.7. Effect of time-step size on the unsteady component lift curves near the end of 
the ramp motion, a/. M c 0 — 0.90, A p = 0.05, a t = 0°, aj — 15°, Re? = 1.52 million. 
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Figure 4.8. Effect of grid refinement on the unsteady component lift curves. = 0.90, 

A p = 0.05, a,- = 0°, a/ = 15°, Re s - 1.52 million. 
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Fig. 4.9 Upper surface streamlines for the steady- 
state canard- wing- body configuration at various 
angles of attack. 
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Figure 4.10. Upper surface instantaneous streamlines at various times for the canard- wing- 
body configuration undergoing unsteady ramp motion. M ^ = 0.90. A p - 0.05, a, = 0°, 
a/ = 12.83°, Re s = 1.52 million. 
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Figure 4.11. Upper surface pressure contours at various times for the canard-wing-bodv 
configuration undergoing unsteady ramp motion. M ^ = 0.90, A p = 0.05, a, = 0 , 

a/ = 12.83°, Re s - 1.52 million. 
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Figure 4.12. Comparison of crossflow total pressure contours between the steady and un 
steady cases at various angles of attack. M ^ = 0.90. A p = 0.05, a, = 0°, a/ - - 12.83° 
Re 5 = 1.52 million. 
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Figure 4.13. Crossflow total pressure contours for the transient flow after the ramp motion 
has ended. A/*, = 0.90, A p = 0.05, a t = 0°, a, = 12.83°, Re- C = 1.52 million. 
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Figure 4.14. Effect of canard on component lift coefficients for two cycles of oscillation. 
Moo = 0.70, k — 0.268, a m — 4.0°, ^4 = 1.0°, Re? — 3.05 million. 
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Figure 4.15. Effect of canard on component pitching moment coefficients for two cycles of 
oscillation. — 0.70, k = 0.268, a m — 4.0°, A = 1.0°, Re c = 3.05 million. 
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Canard-Off 


Canard-On 



Figure 4.16. Effect of canard on lift and pitching moment coefficient curves for third cycle 
of oscillation. M 0 0 — 0.70, k = 0.268, a m = 4.0°, = 1.0°, Rec — 3.05 million. 



Figure 4.17. Lift and pitching moment curves for the canard-off configuration (started from 
steady-state solutions with three cycles of oscillations shown). A/co = 0.70, k = 0.268, 
a m = 4.0°, A — 1.0°, Re d = 3.05 million. 
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Figure 4.18. Comparisons of directly computed and Fourier-coefficient represented lift and 
pitching moment curves. = 0.70, k = 0.268, a m = 4.0°, .4 = 1.0°, Re s = 3.05 

million. 
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Figure 4.19. Comparison of the computed and experimental longitudinal dynamic stability 
parameters for the canard-off configuration. = 0.70, Re s = 3.05 million. 


89 







— A - - Canara-Oi 
• Canard-Off 



a (deg) 

a) Normal force due to pitch rate parameter 



a (deg) 

c) Damping-in-pitch parameter 




Ct (deg) a (deg) 

b) Normal force due to pitch displacement parameter d) Pitch oscillatory longitudinal-stability parameter 


Figure 4.20. Comparison of the longitudinal dynamic stability parameters for the configu- 
ration with and without the canard. M 0 0 = 0.70, Re £ = 3.05 million. 
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Figure 4.21. Effect of the canard on pitching- moment dynamic stability parameters (dif- 
ferences between canard -011 and canard-off values shown). M 0 0 = 0.70, Re c = 3.05 
million. 
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Figure 4.22. Effect of reduced frequency (k) on normal force and pitching moment coeffi- 
dents for the canard-on case. = 0.70, a m = 4.0°, .4 = 1.0°, Re 5 = 3.05 million. 
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Figure 4.23. Effect of reduced frequency (k) on damping-in-pitch and pitch-oscillatory sta- 
bility parameters for the canard-on case. M oo = 0.70, a m - 4.0°, Re 5 = 3.05 million. 
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Chapter 5 

Unsteady Aerodynamics of Moving 
Canards 


Many canard- configured aircraft, including several mentioned in Chapter 1. utilize moving canards 
for longitudinal stability, pitch control and aerodynamic performance. In order to understand the 
influence of moving canards and to demonstrate the current capability to predict the resulting unsteady 
flowfield. time-accurate Navier-Stokes simulations are solved for the configuration undergoing pitch-up 
and oscillatory canard motions. Details of the grid generation, boundary conditions and zonal interfaces 
for the moving canard cases are found in sections 2.4, 2.5, and 2.6. respectively. The freestream flowfield 
conditions used in the canard motion studies are Moc — 0.85 and a Reynolds number based on the 
mean aerodvnamic chord of the wing ( Rgc) of 2.82 million. These conditions match those used for the 
steady-state canard deflection studies (secs. 3.2.2 and 3.3.3), thereby allowing for the direct comparisons 
and observations between fixed and moving canard results. 


5.1 Canard Pitch-Up Ramp Motion 

The deflections of the canard during the pitch-up ramp motions are illustrated in figure 5.1. All canard 
pitch-up cases begin with an initial canard deflection ( 8 Ci ) of 0.0 deg, and final canard deflections (<M 
of 5.0 and 9.25 deg. As with the pitch-up motions of the entire fixed-canard configuration (sec. 4.1). 
the canard motions begin and end instantaneously. Various canard pitch rates ( A pc ) are simulated 
in this section. The time histories of force and moment are again given in non-dimensional time 
(t raJ = y'otf that during the portion when the canard is moving, 8 C = t. 

5.1.1 Unsteady lift and pitching moment 

The time histories of lift and pitching moment coefficients are given in figure 5.2 for the canard con- 
figuration at a = 4.27 deg. Canard deflection angles as a function of non-dimensional time (t) are 
also indicated in the figure. Starting from a steady-state condition at 8 C = 0.0 deg, the canard begins 
to pitch up at time, t = 0 deg with a pitch rate of Ap C = 0.05. The lift curve for o 4.2 1 deg 
(fig. 5.2(a)) shoves that an increase in total lift is purely transitory during the unsteady aerodynamic 
portion of the canard motion. In fact, the initial total lift at t = 0 deg (<5 C = 0.0 deg) and the 
final total lift for t > 8 deg (S c = 5.0 deg) are basically unchanged. This result is consistent with the 
steady-state computations at various fixed canard deflections shown earlier (fig. 3.14), which illustrated 
an increase in canard lift and a decrease in wing lift as a function of increasing canard deflection angles. 


95 



Figure 5.2(b) shows the effects of the canard motion on the pitching moment coefficients for the 
configuration. Steady-state results of section 3.2 found that the wing portion pitching moment is 
relatively insensitive to the various fixed canard deflection angles. The results in these two figures 
indicate that this effect is extended to unsteady canard ramp motions as well. For the configuration in 
this study, the pitching moments for the entire configuration at a = 4.27 deg (fig. 5.2(b)) is purely 
dependent on the canard portion contributions. 

Although the pitching moments of the configuration undergoing canard pitch-up motion are very 
similar to the corresponding steady-state characteristics, the unsteady nature of the lift curve is of 
particular interest. Figure 5.2(a) also shows that although the canard portion lift responds almost 
immediately to the canard pitch-up motion, a significant delay is observed before the wing portion lift 
responds to the effects of the canard motion. This delay is due to the time required to convect the 
effects of the canard downstream to the wing. The length of the delay observed in the figure, which is 
on the order of A t « 3 deg (or 0.05 rad.), correlates well with the pitch rate of A pc = 0.05 for this 
angle-of-attack case. In essence, for A pc = 0.05, At — 0.05 rad. is the length of time for a freestream 
fluid particle to travel the distance c, which also represents the approximate average distance between 
the canard and the wing. The transient increase in total lift is primarily due to this delay in the wing’s 
response to the increase in canard lift. 

The transient nature of the configuration lift curve and, specifically, the delay observed in the wing 
response is a direct function of the canard pitch rate. Figures 5.3 and 5.4 show' the effects of canard 
pitch rate, A pc , on lift and moment coefficients, respectively, for the configuration at a = 4.27 deg. 
As pitch rate is increased, both the magnitudes of the canard lift and the non-dimensional time delay 
in the wing response increase (figs. 5.3(b) and (c)). The combination of these two effects produces a 
higher peak lift coefficient for the entire configuration (fig. 5.3(a)). The results of figure 5.4 illustrate 
that canard pitch-rate effects on pitching moment are considerably less significant than for lift. 

For the higher pitch rate cases given in figures 5.3 and 5.4, a constant physical time-step size results 
in a higher non-dimensional time step. However, refining the step sizes can be used to demonstrate 
temporal convergence. Figure 5.5 show's the effect of time-step size on the lift and pitching moment 
predictions of the computations performed for A pc = 0.20. The excellent comparison between the 
computational simulations at various step sizes, particularly during the transient period after the canard 
motion has stopped (t > 5 deg), confirms the adequate time-step resolution for the canard pitch-up 
ramp cases. 

5.1.2 Canard-wing vortex interaction 

The results given in the preceding section can be further understood by examining the canard-wing 
vortex interaction during the canard pitch-up maneuver. For easier visualization of stronger vor- 
tices, a canard pitch-up simulation is performed for the configuration at the higher angle of attack 
of a = 8.0 deg. Figure 5.6 illustrates the time histories of lift and pitching moment coefficients 

for this case. Similar to the lower angle-of-attack case given earlier in figure 5.2, the current case at 
a = 8.0 deg shows a transitory increase in total lift which is due to the delay in wfing lift response. 
Figure 5.6(b) also show's that pitching moment characteristics during the canard motion are similar to 
the lower angle-of-attack case. 

Comparisons of total pressure contours at constant streamwise locations are given in figure 5.7 
for the instant wffien the canard is deflected at S c = 5.0 deg during the canard pitch-up motion 

( A pc = 0.05). The steady- state results of the static canard at the corresponding freestream condition 
and canard deflection angles are also shown in the figure. The crossflow' total pressure contours illustrate 
a slightly weaker canard vortex formation for the unsteady cases (xjc — 2.0) in comparison to the 
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steady cases. The weaker canard vortex at a given deflection angle is primarily due to the induced 
velocities at the canard leading edge from the canard rotation about its pitch axis. Further downstream 
at the same instant in time, the weaker canard vortex and the effects from convection of the canard's 
influences produce a less developed canard vortex over the wing surface ( xjc = 2.8 and 3.5). After 
the canard motions stops, the effects of the canard vortex on the wing flowfield will approach the 
steady-state results. 

To better visualize the effects of the wing flowfield development, the time scales and canard effects 
are increased by starting with the lower angle of attack case (ct = 4.27 deg) and continuing the 

canard pitch-up motion to a higher deflection angle. Figure 5.8 shows the lift and pitching moment 
curves for the configuration undergoing a canard pitch-up from <5 Ci = 0 deg to 6 CJ = 9.25 deg with 
a non-dimensional pitch rate of A pc = 0.05. The increase in total lift and the delay in wing response 
are both clearly evident in figure 5.8(a). For the canard pitching to 8 C = 9.25 deg. the downwash of 
the canard causes the wing to generate very little lift at the positive a of 4.27 deg. This effect was 
also observed in the deflected static canard cases of section 3.2.2, which showed steady-state vortices 
on both the wing upper and lower surfaces. 

Traces representing instantaneous streamlines at seven different times are given in figure 5.9 for the 
8 C{ — 9.25 deg canard pitch-up case. As illustrated in the legend of figure 5.9, the first four instances 
in time (A-D) are during the canard motion (6 C < <5 C/ ), and the last three instances (E-G) are during 
the transient phase as the flowfield approaches a steady-state condition. Visualization of the canard 
and wing vortices is performed by starting the traces near the upper- and lower-surface leading edges 
of the wing, and the upper-surface leading edge of the canard. 

As the canard deflection begins to increase, minimal initial influence of the canard on the wing 
vortex structure is shown (time positions A and B). As the canard deflection approaches 8c ~ 7 deg, 
time position C shows the initial formation of a leading edge vortex on the wing lower surface. Just 
before the canard motion stops at time position D, the wing lower surface vortex becomes more or- 
ganized. and the canard leading-edge vortex shows evidence of vortex breakdown. Time positions 
E-G show the completion of the vortex structure development while the canard deflection is fixed at 
S c = 8 Cf = 9.25 deg. The formations of the lower surface vortex on the inboard of the wing and the 
upper surface vortex on the outboard wing portion are clearly evident. 

Crossflow visualizations of total pressure contours are given in figure 5.10 for the canard and wing 
regions. Figures 5.10(a)-(c) correspond to time instances A-C from figures 5.9(a)-(c) during the canard 
ramp motion. In figure 5.10, the development of the canard vortex and its convection over the wing 
surface is observed. Furthermore, the initial formation of the wing lower surface vortex is noted in 
figure 5.10(c) (wing region) and correlates with the instantaneous streamlines of figure 5.9(c). 

In order to more clearly visualize the development of the wing vortex structure, crossflow total 
pressure contours are given in figure 5.11 at the same xjc — 3.3 for five time positions starting with 
just prior to when the canard deflection angle has reached 8 Cf . As time increases, the development of 
the canard vortex over the wing and the wing leading-edge vortices on the upper and lower surfaces 
are observed. Since the contour levels of total pressure are evenly distributed, the relative strengths 
of these vortices are indicated in figure 5.11 as an approximate loss in total pressure at the vortex 
cores. Although at slightly different conditions, a qualitative comparison can also be made with the 
static canard deflection cases given in figure 3.26, w-here the canard, wing upper surface and wing lower 
surface vortices wore also observed in the wing flowfield. As illustrated earlier in figures 5.2 and o.6 of 
section 5.1.1, and currently in figures 5.9 and 5.11, the time required for the full development of the 
wing flowfield and vortex structure allows for the transient increases observed in the total lift of the 
configuration. 
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5.2 Canard Pitch Oscillation 


The investigation of the unsteady aerodynamics of canard pitch oscillations can lead to an evaluation 
of the dynamic stability' characteristics of the moving canard control surface. Furthermore, previous 
studies (refs. 26-28) have experimentally examined the potential of using canard oscillations to promote 
wing performance in extreme maneuver conditions. Specifically, these studies explored the potential of 
using high frequency canard oscillations as an active control method of delaying wing stall. Although 
unsteady experiments may be able to evaluate the effectiveness of such techniques, computational 
simulations are required to fully understand the physics of such unsteady flow phenomenon and. more 
importantly, to optimize the use of the oscillating canard. This section will present the results of canard 
pitch oscillations at the pre-stall conditions of interest in this study. 

Canard pitch oscillation simulations are performed at a pitch amplitude [A c ) of 2 deg, and reduced 
frequencies (k c ) of 1.0. Mean angles of attack for the configuration of 4.27 deg and 8.00 deg are 
presented. As with the rigid-body oscillation studies of sections 4.2 and 4.3. high values of reduced 
frequencies are simulated to minimize computational expense. 

5.2.1 Unsteady lift and pitching moment 

The effects of the oscillating canard on the configuration lift and pitching moment characteristics 
are given in figures 5.12(a) and (b), respectively', as a function of non-dimensional time {ui c t) for the 
configuration at. a = 4.27 deg. The canard deflection angles with respect to ui c t are also indicated 
in the figure. From figure 5.12. it is evident that the time histories of the canard lift and moment are 
closely in-phase with the canard oscillation. On the other hand, the wing contributions to both lift 
and moment exhibit significant shifts in phase and appear to be the primary contributor to the phase 
differences between the canard motion and the total lifts and moments. These effects of the oscillating 
canard on the canard and wing contributions are consistent with the earlier findings of the canard 
pitch-up ramp motion studies (sec. 5.1). which also found almost instantaneous canard response to 
changes in deflection angle 6 C . 

Figure 5.13 shows that a periodic time response is achieved within two cycles of canard oscillations in 
spite of the fact that the wing responses are due to the canard motions rather than direct wing motions. 
Note also that, the significant contributions of the wing responses to the overall phase differences in 
both lift, and moment are observed in figures 5.13(a) and (b), respectively. In fact, figure 5.13(b) shows 
that the phase differences between the overall pitching moment and the canard deflections are almost 
entirely due to the wing response. Furthermore, it is significant, to observe that the relatively small 
amplitude oscillations of the canard produce relatively large wing responses in pitching moment. The 
amplitude of wing responses due to the pitch oscillations are comparable to those resulting from the 
much larger canard deflections of the canard ramp motion cases (fig. 5.8). Correspondingly, the wing- 
lift response due to the canard oscillations in figure 5.13(a) also exhibits a relatively high amplitude 
for the small canard deflection angles. These results indicate that the high frequency nature of the 
canard oscillations has a significant influence on the wing response characteristics, even at relatively 
small pitch amplitudes. 

Similar trends are observed for the lift and moment characteristics of the configuration at 
a = 8.0 deg. Figures 5.14 and 5.15 show that the canard oscillation produces nearly in-phase re- 
sponses for the canard region and significant phase shifts for the wing responses in both lift and pitching 
moment. Amplitudes of the total configuration as well as the canard and wing region responses are 
comparable to the lower angle-of-attack case. 
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5.2.2 Canard and wing vortex structures 

Visualizations of the canard and wing vortex structures are given in figures 0.16 and 5.1 1 , respectively, 
for the canard oscillation case of a = 4.27 deg. For both figures, crossflow visualization of normalized 
total pressure contours are shown at different u) c t values, as indicated in the legend at the end of each 
figure. Figure 5.16 shows that the strength and location of the canard leading-edge vortex changes 
throughout the cycle. Early in the cycle (a-c), the canard vortex appears weaker and is increasing in 
strength. Near the middle to the end of the cycle (d-i), the canard vortex appears to move inward 
towards the root of the canard' and then begins to weaken again. The effect on the wing of these 
changes in the canard vortex characteristics is shown in figure 5.17. Over the stationary wing surface, 
the converted canard vortex changes in strength, location and structure throughout the oscillation 
cycle. It is noted from figure 5.17 that the wing vortex strength and location appear to be relatively 
unchanged. However, further analysis of the wing surface pressures shows that the characteristics of 
the wing vortex are indeed affected by the oscillating canard. 

Figure 5.18 illustrates the effect of the oscillating canard on the wing surface pressures at the 
55 percent and 75 percent semispan stations. Each streamwise pressure distribution curve represents 
the instantaneous pressure on the wing at the same corresponding time instance (A, C. E. G. and I) 
indicated earlier in figures 5.16(a)-5.17(i). On the inboard portion of the wing, the canard oscillation 
clearly affects the formation of the wing leading-edge vortex. The effect on the lower-surface pressure 
distribution near the leading edge indicates that the local effective angle of attack of the wing is changing 
as well. Further outboard, figure 5.18 shows that there are significant changes in the upper-surface 
pressure peaks which are generated by the w r ing vortex. Starting from the beginning of the cy cle, the 
suction peak weakens and appears to move aft (from A— C). Later in the same cycle, the peak moves 
forward towards the leading edge as suction again increases near the end of the cycle (E-I). Although 
the wing upper surface at 75 percent semispan shows considerable change, the local effective angle 
of attack at this station is observed to be relatively constant from the nearly constant lower surface 
leading-edge pressures. 

The effects of the canard oscillation on vortex structures are given in figures 5.19 and 5.20 for the 
configuration at Ci — 8.0 deg. At these conditions, the amplitude of the canard oscillations is a 

smaller fraction of the mean angle of attack. Therefore, over the canard region (fig. 5.19), the changes 
in the canard vortex characteristics appear to be more subtle than for the earlier a = 4.27 deg case. 
How’ever, results over the wing in figure 5.20 show that the converted canard vortex does indeed change 
in strength and structure throughout the canard oscillation cycle. Similar to the earlier a = 4.27 deg 
case, figure 5.21 shows that the wing vortex at a = 8.0 deg, as characterized through surface pressures, 
is also affected by the canard oscillations. The local effective angle of attack of the wing is again most 
influenced in the inboard portion of the wing. 
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a) Canard ramp to 5.00 deg. 



b) Canard ramp to 9.25 deg. 


Figure 5.1. Schematic of canard pitch-up ramp motion from 6 Ct = 0° to 6 Cf = 5° and 9.25°. 
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a) Lift coefficient 



b) Pitching moment coefficient 


Figure 5.2. Time history of unsteady lift and pitching-moment responses due to the canard 
pitch-up ramp motion for the configuration at a = 4.27°. M c 0 = 0-85, A pc = 0.05, = 0°, 

6 C/ = 5°, Re 5 = 2.82 million. 
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Figure 5.3. Effect of canard pitch rate on the time history of unsteady lift responses due to 
the canard pitch-up ramp motion. M 0 0 = 0.85, a — 4.27°. S Ci ~ 0°, 6 CJ = 5°, Re d = 2.82 
million. 
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Figure 5.4. Effect of canard pitch rate on the time history of unsteady pitching moment 
responses due to the canard pitch-up ramp motion. M 0 0 = 0.85, a = 4.27°, 6 Ct = 0°, 

6 Cf — 5°, Re 5 = 2.82 million. 
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Figure 5.5. Effect of time-step size on unsteady lift and moment responses due to the 
canard pitch-up ramp motion for A pc = 0.20. M = 0.85, a = 4.27°, <5 C| = 0°, <5 C/ = 5°, 
Re 5 = 2.82 million. 
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a) Lift Coefficient 



b) Pitching Moment Coefficient 


Figure 5.6. Time history of unsteady lift and pitching-moment responses due to the canard 
pitch-up ramp motion for the configuration at a = 8.0°. M 0 0 = 0.85, A pc — 0.05, = 0 , 

6 Cf = 5°, Re 5 = 2.82 million. 
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STEADY-STATE 
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x/c = 2.8 




x/c = 3.5 

Figure 5.7. Comparison of crossflow total pressure contours between the steady and un- 
steady cases at various streamwise locations for the configuration at a = 8.0° and S c = 5.0° 
(unsteady case given at the instantaneous 6 C for t — 5.0°). = 0.85, A pc = 0.05, 

S Ci = 0°, S CJ — 5°, Rti — 2.82 million. 
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a) Lilt Coefficient 



Figure 5.8. Time history of unsteady lift and pitching-moment responses due to the canard 
pitch-up ramp motion from 6 Ct = 0° to 6 Cf — 9.25°. M oo — 0.85. a = 4.27°, .T pc = 0.05, 
Re 5 = 2.82 million. 
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UPPER SURFACE 


LOWER SURFACE 



a) Time Instance A: t = 2.33 deg; 5 C = 2.33 deg 



b) Time Instance B: t = 4.66 deg; 8 C = 4.66 deg 



c) Time Instance C: t = 7.00 deg; 5c = 7.00 deg 



d) Time Instance D: t = 9.19 deg; 5c = 9.19 deg 

Figure 5.9. (Part-1) Upper and lower surface instantaneous streamlines for the canard 
pitch-up ramp motion from 6 Ci - 0° to 6 Cf = 9.25°. M 0 0 — 0.85, a = 4.27°, A pc — 0.05, 
Re? — 2.82 million. 
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e) Time Instance E: t = 12.00 deg; 5c - 9.25 deg 



f) Time Instance F: t = 13.67 deg; 5c - 9.25 deg 



g) Time Instance G: t = 18.33 deg; 5c - 9.25 deg 
LEGEND 



iim» (4ag) 


Figure 5.9. (Part-2) Upper and lower surface instantaneous 
pitch-up ramp motion from 6 Cx - 0° to 6 Cf = 9.25 . A/oo - 


streamlines for the canard 
0.85, a = 4.27°, A pc = 0.05, 


Re s = 2.82 million. 
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CANARD REGION (x/c = 2.0) 


WING REGION (x/c = 3.3) 





c) t = 7.00 deg 


Figure 5.10. Crossflow total pressure contours over the canard (x/c = 2.0) and wing ( x/c = 
3.3) during the canard pitch-up ramp motion from 6 Ci = 0° to S Cf = 9.25°. = 0.85, 

a = 4.27°, A pc = 0.05, Re e = 2.82 million. 
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c) t = 12.00 deg (Canard at 9.25 deg) 


Figure 5.11. Crossflow total pressure contours over the wing (x/c = 3.3) after the canard 
pitch-up ramp motion from 6 Ct = 0° to 6 Cf = 9.25°. A/,*, = 0.85, a = 4.27°, A pc — 0.05, 
Re s - 2.82 million. 
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a) Lift Coefficient 



b) Pitching Moment Coefficient 


Figure 5.12. Time history of unsteady lift and pitching-moment responses due to the canard 
pitch oscillation motion for the configuration at a = 4.27°. M ^ = 0.85, A c — 2.0°, 

k c — 1.0, Re d = 2.82 million. 


112 





a) Lift Coefficient 



b) Pitching Moment Coefficient 


Figure 5.13. Lift and pitching-moment curves for the configuration undergoing canard pitch 
oscillation motion at a = 4.27°. M 0 0 = 0.85, A c = 2.0°, k e = 1.0, Re- C = 2.82 milUon. 
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a) Lift Coefficient 
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Figure 5,14. Time history of unsteady lift and pitching-moment responses due to the canard 
pitch oscillation motion for the configuration at a = 8,0°. = 0.85, A c — 2.0°, k c = 1.0, 

Re d = 2.82 million. 
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Figure 5.15. Lift and pitching-moment curves for the configuration undergoing canard pitch 
oscillation motion at a — 8.0°. Mqq — 0.85, A c = 2.0°, k c = 1.0, Res — 2.82 million. 
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c) co c t = 90.0 deg. *) “c 1 = 223 - 2 de g- 


Figure 5.16. (Part-1) Crossflow total pressure contours over the canard (x/c - 2.0) for the 
canard pitch oscillation motion at a = 4.27°. M 0 0 = 0.85, A c = 2.0°, k c = 1.0, Re e = 2.82 
million. 
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g) co c t = 270.0 deg. 
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h) oo c t = 313.2 deg. 



i) co c t = 360.0 deg. 


Figure 5.16. (Part-2) Crossflow total pressure contours over the canard (x/c = 2.0) for the 
canard pitch oscillation motion at a = 4.27°. M <*> 0.85, /l c - ‘2-0°, = CO, Re s = 2.82 

million. 




c) co c t = 90.0 deg. f) co c t = 223.2 deg. 


Figure 5.17. (Part-1) Crossflow total pressure contours over the wing ( x/c = 3.5) for the 
canard pitch oscillation motion at a = 4.27°. = 0.85. A c - 2.0°, k c = 1.0, ifec = 2.82 

million. 
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LEGEND 



h)oo c t = 313.2 



i) <o c t = 360.0 deg. 

Figure 5.17. (Part-2) Crossflow total pressure contours over the wing (x/c = 3.5) for the 
canard pitch oscillation motion at a = 4.27°. M 0 0 = 0.85, A c = 2.0°, k c — 1.0, Re £ — 2.82 












a) co c t = 0 deg. 


d)co c t= 133.2 deg. 



b) co c t = 43.2 deg. 


e) G) c t = 180.0 deg. 



C)0) c t = 90.0 deg. 0co c t = 223.2 deg. 

Figure 5.19. (Part-1) Crossflow total pressure contours over the canard (x/c = 2.0) for the 
canard pitch oscillation motion at a = 8.0°. Moc = 0.85, A c — 2.0°, k c = L.0, Re^ = 2.82 
million. 
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LEGEND 



h) to c t = 313.2 deg. 



i) co c t = 360.0 deg. 


Figure 5.19. (Part-2) Crossflow total pressure contours over the canard ( x/c = 2.0) for the 
canard pitch oscillation motion at a — 8.0°. = 0.85, A c = 2.0°, k c = 1.0, Re s — 2.82 

million. 
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c)co c t = 90.0 deg. f)< 0 c t = 223.2 deg. 

Figure 5.20. (Part-1) Crossflow total pressure contours over the wing ( xjc — 3.5) for the 
canard pitch oscillation motion at a — 8.0°. — 0.85, A c = 2.0°, k c = 1.0, Rtc — 2.82 

million. 
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i) <o c t = 360.0 deg. 

Figure 5.20. (Part-2) Crossflow total pressure contours over the wing (x/c = 3.5) for the 
canard pitch oscillation motion at a = 8.0°. Moo = 0.85, A c = 2.0°, k c = 1.0, Re 5 = 2.82 



a) 55% Span Station 



LEGEND 



CD. t (deg) 

Figure 5.21. Effect of canard pitch oscillation motion on wing surface pressures at a 
Mm = 0.85, A c = 2.0°. k c = 1.0, Re 5 = 2.82 million. 






Chapter 6 


Conclusions 


A detailed numerical investigation has been performed for the steady and unsteady flowfield about a 
static and dynamic canard- wing-body configuration. Computational technology has been developed 
and enhanced to perform a series of time-accurate Navier-Stokes simulations about the close-coupled 
canard configuration. The computational methods presented in this study are capable of accurately 
predicting the unsteady viscous flowfield for the canard-wing-bodv configuration with various fixed 
canard positions and with moving canards. 

6.1 Validation of the Computational Method 

The accuracy of the simulations is demonstrated by extensive comparisons with available steady and 
unsteady’ experimental data. Favorable comparisons with steady-state experimental surface pressures, 
forces, and moments for the configuration with various fixed canard positions verify that the present 
method is capable of accurately predicting the steady-state flowfield about such configurations. For 
low angles of attack (up to approximately 6 deg), a relatively’ coarse (baseline) grid was found to be 
effective in accurately computing surface pressures, forces and moments. However, a more refined 
grid was shown to be critical for higher angle-of-attack cases where boundary -layer separation, strong 
shocks, and vortex breakdown dominate the flow characteristics. 

Comparisons of experimentally measured dynamic stability parameters with computations of the 
configuration undergoing pitch oscillations also demonstrate accurate predictions of the unsteady 
canard- wing- body flowfield. An analysis of the effect of reduced frequency was also performed and 
demonstrates the ability to accurately predict dynamic stability with considerably higher computa- 
tional reduced frequencies. As a result, significant reductions in required computational resources 
allows the use of state-of-the-art computational fluid dynamics (CFD) for the prediction of dynamic 
stability about complex canard configurations. For those cases where experimental data was unavail- 
able, a spatial and time-step refinement study was conducted to verify the convergence of the unsteady 
computations. 

6.2 Canard-Wing-Body Steady Flowfield 

For the steady-state mid-canard flowfield, the canard-induced downwash was found to weaken or delay- 
formation of the wing leading-edge vortex. At a « 4 deg, the formation of the wing leading-edge vortex 
is delayed to the canard-tip span-line of the wing. At a « 8 deg and 12 deg. a weakened leading-edge 
vortex was observed on the inboard portion of the wing. The wing vortex on the outboard span was 
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shown to be influenced by both the canard vortex and the aforementioned canard downwash. The 
presence of the canard was also shown to eliminate the wing vortex breakdown which was evident in 
the canard-off case at a ~ 12 deg. These results confirm and provide an explanation for the canard’s 
potential for delaying wing vortex breakdown which has been observed in experimental studies. 

At low to moderate angles of attack, the fixed high-canard configuration exhibits improved steady- 
state lift and drag over the mid-canard case. However, due to the interaction between the low-canard 
vortex and the wing surface, unfavorable lift and drag characteristics for the low-canard configuration 
are evident as angle of attack is increased. At low angles of attack, the low-canard vortex passes under 
the wing surface and can induce lower pressures on the wing lower surface. At higher angles of attack, 
the canard vortex is split by the wing surface into upper and lower vortices. Toward the outboard 
wing, the wing vortex induces a large inward movement of the canard upper-vortex. 

When the fixed canard is positively deflected, a pronounced effect on aerodynamic performance 
parameters is observed, particularly for wing lift and configuration pitching moment. Visualization of 
the canard-wing-body flowfield shows a complex flow structure consisting of several interacting canard 
and wing vortices. Specifically, at significantly large canard deflection angles, effective wing angle of 
attack transitions from negative (inboard) to positive (outboard) values causing the formation of both 
lower and upper surface leading-edge vortices. Furthermore, due to the strength of these vortices, 
secondary vortices are present indicating strong viscous effects. 

These findings both correlate with, and provide an explanation for, much of the experimental 
results presented in the literature. These results also suggest that detailed investigations are required 
to determine optimum canard positioning for given configurations. The accurate prediction of such 
steady-state effects are needed in the design and optimization of future close-coupled canard aircraft. 


6.3 Canard-Wing-Body Unsteady Flowfield 

To understand the predicted unsteady flowfield, an analysis of the unsteady canard-wing-body aero- 
dynamics with emphasis on the canard-wing vortex interaction was performed. Dynamic conditions in 
this study included rigid-body pitch-up and pitch oscillation motions of both the entire configuration 
and the independent canard. 

For the entire configuration undergoing pitch-up ramp motions, computations with and without the 
canard showed increased dynamic lift and improved lift-to-drag ratios for the configuration with the 
canard. Results also indicated the strong influence of the relative positioning of the canard and wing 
to the pitch axis. For the canard located forward of the pitch axis, the ramp motion at high pitch rates 
decreases the effective angle of attack of the canard and delays the formation of the canard vortex. 
However, as angle of attack increases and the ramp motion ends, the canard vortex and its interaction 
with the wing develops rapidly. As time increases, the configuration flowfield approaches a relatively 
steady solution. 

During pitch oscillations, the fixed canard was found to have a pronounced effect on the unsteady 
aerodynamics associated with the wing-body configuration. Significant changes in amplitude and shifts 
in phase of the lifts and moments were observed with the addition of the canard. Using the results 
of the numerical pitch oscillation simulation, an analysis of the time histories of normal force and 
pitching moment coefficients led to the prediction of longitudinal dynamic stability parameters. These 
parameters were utilized to evaluate the effects of the canard on dynamic stability. The current study 
found that the presence of the canard increases the pitch damping of the configuration and decreases 
the oscillatory stability. 
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The current computational method has been further extended to investigate the unsteady aero- 
dynamic effects of moving canards. Both canard pitch-up ramp and pitch oscillation motions were 
simulated using the developed dynamic gridding, zonal interfacing and boundary condition capabili- 
ties. The canard pitch-up ramp motion simulation showed that a transient increase in overall lift was 
the direct result of the time delay (or lag) in wing response from the convection of the canard ef- 
fects. Furthermore, the vortex analyses illustrated the unsteady nature of the canard and wing vortex 
structures, particularly with respect to the formation of the upper- and lower-surface wing vortices. 

The investigation of canard pitch oscillations demonstrated that high frequency, small amplitude 
oscillations can produce significant lift and moment responses of the wing. It was also observed that, 
similar to the canard ramp motion results, the canard responses were nearly in-phase with the canard 
motion while the wing responses showed significant phase shifts. Although the vortex visualizations 
indicated minimal changes in wing vortex characteristics, analyses of the wing surface pressures showed 
significant changes in the local angle of attack on the inboard wing and overall changes in wing vortex 
influences. 

6.4 Recommendations for Future Studies 

The current study has characterized the complex steady and unsteady aerodynamic flow-field of a canard 
configuration, including effects of both fixed canards at various positions and independently mo\ing 
canards. Future studies should focus on utilizing the tools and findings which have been presented 
to optimize the use of canards for both cruise and maneuver aerodynamic performance. Such studies 
could include, but are not limited to: 

1. Optimization of fixed canard positions (vertical, longitudinal, and deflection) for trimmed cruise 
performance; 

2. Development of control-law based canard motions for enhanced maneuver performance; 

3. Investigation of non-symmetric canard positions or motions for lateral stability and control; 

4. Investigation of high-frequency canard oscillations for improved wung performance in pre- and 
post-stall environments; 

5. Evaluation of complete static and dynamic stability characteristics of a canard configuration for 

stability and controls applications; and 

6. Coupled controls- fluids simulation of a dynamic maneuver condition. 
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Appendix A 
Governing Equations 


(Note: For further details of the governing equations , numerical algorithm, and turbulence modeling, 
the reader is referred to ref. f'3 and the appropriate original publications which are cited.) 

The governing equations solved in this study are nondimensionalized by free-stream quantities and 
are transformed to the computational domain (£. r},() so as to preserve the strong conservation form 
of the equations. The strong conservation law form of the Navier- Stokes equations are used for shock 
capturing purposes. The equations in Cartesian coordinates in non-dimensional form can be written 


as 


where 
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The Cartesian velocity components u. v, and w are nondimensionalized by a x (the free-stream 
speed of sound), density p is nondimensionalized by p x ; and the total energy per unit volume e is 
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(A -4) 


nondimensionalized by poetic - Pressure can be found from the ideal gas law as 

p =( 7 - l)[e - 0.5p(u 2 + v 2 + u- 2 )] 

and throughout 7 is the ratio of the specific heats. Also, n is the coefficient of thermal conductivity, 
p is the dynamic viscosity, and A from the Stokes' hypothesis is - 2 / 3 /; . The Reynolds number is Re 
and the Prandtl number is Pr. 

To enhance numerical accuracy and efficiency and to handle boundary conditions more easily, the 
governing equations are transformed from the Cartesian coordinates to general curvilinear coordinates 
where 


T — t 

£ = £{x,y,z,t) (A - 5) 
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C = C (x,y,z,t) 


The resulting transformed equations are not much more complicated than the original Cartesian 
set and can be written in nondimensional form as 
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\\ = Q + Cx u + C y v + Q w 

where £7, K and W are contravariant velocity components written without metric normalization. The 
viscous flux terms are given by 
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where the components of the shear-stress tensor and heat-flux vector in nondimensional form were 
given in equation A-3. Here, the Cartesian derivatives are expanded in £, rj, £ space via chain-rule 
relations such as 


Ux — £x^£ "F Vx^rj ~F Cx^£ 


Finally, the metric terms are obtained from chain-rule expansion of x^. y 7? , etc., and solved for £y, 


etc., to give 
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and 
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In high Revnolds number flows, the viscous effects are confined to a thin layer near rigid boundaries. 
In most cases, there are only enough grid points to resolve the gradients normal to the body by clustering 
the grid in the normal direction, and resolution along the body is similar to what is needed in inviscid 
flow. As a result,, even though the full derivatives are retained in the equations, the gradients along the 
body are not resolved unless the streamwise and circumferential grid spacings are sufficiently small. 
Hence, for many Navier-Stokes computations, the viscous derivatives along the body are dropped. This 
leads to the thin-layer Navier-Stokes equations. 

The thin-layer model requires a boundary layer type coordinate system. For example, if £ and rj 
directions are along the body and the viscous derivatives associated with these directions are dropped, 
the terms in C only are retained and the body surface is mapped onto a constant < surface. As a 
result, equation A-6 simplifies to 

d r Q *F d?E -f- djfF + OqG — Re *<9^5 (A — 11) 
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To generalize the coordinates of complex geometries, the thin-layer approximation can be made in 
all three coordinate directions and. as a result, only cross-derivative viscous terms are neglected. It 
should also be emphasized that the thin-laver approximation is valid only for high Reynolds number 
flows. Very large turbulent eddy viscosities would invalidate the thin-layer assumption. 
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As mentioned in section 2.1. the algorithm used to solve the Reynolds-averaged Navier-Stokes equa- 
tions is due to Beam- Warming. The basic Beam- Warming algorithm is first- or second-order accurate 
in time and second- or fourth-order accurate in space. The equations are approximately factored (spa- 
tially split) which, for a given time iteration, reduces the process to three one-dimensional problems. 
Due to the second-order central-difference operators employed, the algorithm produces block tridiag- 
onal systems for each spatial dimension. The stability and accuracy of the numerical algorithm is 
described bv Beam and Warming (ref. 40). According to the linear analysis, the numerical scheme is 
unconditionally stable in two dimensions but in actual practice time step limits are encountered because 
of the nonlinear nature of the equations. However, this limitation is much less stringent than compa- 
rable explicit schemes. In three dimensions the algorithm is unconditionally unstable, but stability is 
maintained by the addition of artificial dissipation terms. 

The finite-difference algorithm due to Beam and W’arming applied to equation A- 11 results in the 
following approximate factorization: 

(7 + %i n + Df ] ){I + h8 n B n + D f>) 

x(J + h6(C n - hRe~ 1 6 c J- 1 M n J+ D[ 2) ) AQ n = R n (B - 1) 

= -Ar(6cE n + 8r,F n + 8 c G n - Re~ l b Q S n ) - D w Q n 

where 8 is the central-difference operator and A and V are forward and backward-difference operators, 
e.g.. 

b'Q = [Q(£ + A£, 77 . 0 - Q{£, ~ A£, 77 , C)]/ 2A £ 

A <Q = [Q(t + *S,V,0-Q&rhQ]/Mi ( B -2) 

VcQ = 

Indices denoting spatial location are suppressed and h = At corresponds to first-order time- 
accurate Euler Implicit and h = Ar/2 to second-order time-accurate Trapezoidal Rule. D/ .Dr," . 

and DV 2 ' 1 are the implicit and D, <n is the explicit smoothing operator which are given in detail in 

S 

reference 43. A 

The Jacobian matrices A". B n . and <7* are obtained by linearizing the flux vectors E n . F n , and G n 

in time such that 
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where indices denoting spatial location are suppressed again and 


A-°*. B= d -Z . 

dQ dQ 


C = 


do 

dQ 


(B — 4) 


are the flux Jacobian matrices. These flux Jacobians and the viscous coefficient matrix M. which comes 
from the time linearization of the viscous vector 5 n+1 , are documented in reference 43. 

Block tridiagonal-matrix inversions constitute the major portion of numerical work associated with 
the standard Beam- Warming algorithm. Equations A-6 are a coupled set of five equations and thereby 
produce a (5 x 5) block-tridiagonal structure for the implicit operators of equations B-l. The diagonal 
version of the standard algorithm due to Pulliam and Chaussee (ref. 41) overcomes this difficulty. In 
this algorithm, rather than inverting block-tridiagonal matrices in each direction, scalar pentadiagonal 
matrices are inverted. This is computationally more efficient. 

The Jacobian matrices. .4, B, and C have a set of eigenvalues and a complete set of distinct 
eigenvectors. Similarity transformations can be used to diagonalize A, B. and C. 

A = Tc\(rr\ b = T ri A ri T r ;\ c = r c A c r - 1 (B- 5 ) 

where 
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where c is the speed of sound (c~ - 7 p/p), and for example Ac reads as 
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The similarity transformation matrices Tc, T n , and their inverse matrices are given in reference 43. 
Relations exist between Tc. T v , and Tq of the form 

N = Tr 1 ^. N- 1 = T~%, P - T~ l T c . P~ ] = (B - 8) 


where 
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After applying the similarity transformations of equation B-5 and identities given in equation B-8 into 
equation B-l and exchanging the smoothing operators with new ones, the diagonal form of the standard 
algorithm reads 

2>(J + h6s Ac - hDi\z)N(I + h8 v A v - hD l \ r ,)P(I + - hD l \^)T~ 1 /AQ n = R n (B - 10) 

The spatial accuracy of the standard and diagonalized algorithms for steady-state problems (i.e., 
/\Q n — * o as n — ► oc) is determined by the type of differencing in forming R n . Since the modification 
that produces the diagonal algorithm does not effect R n , both schemes will have the same steady- 
state solution assuming that the steady-state solution is independent of the convergence path; i.e., that 
the steady state is unique. For unsteady simulations, care must be exercised in ensuring small enough 
time-step sizes for accuracy. When the diagonalized algorithm is used, the author strongly recommends 
time-step refinement studies to evaluate temporal accuracy of the unsteady simulations. 

For constant coefficient matrices A, B, and C 1 , the diagonal algorithm reduces to the standard 
algorithm because the eigenvector matrices are also constant. Therefore, the linear stability analysis 
of Beam and Warming also holds for the diagonal algorithm. 
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The apparent turbulent stresses in compressible flow can be written as 

(Tij)turb = PU-^Uj 


(C-l) 


and the apparent turbulent heat flux components as 

Q 

~(V • q )turb = -far ( pC P T ' U ^ 


(C-2) 


In order to predict turbulent flows by applying finite-difference methods to the Reynolds equations, 
it is necessary to make assumptions for the terms in equations C-l and C-2. Boussinesq (18/ / ) suggested 
that the apparent shearing stresses might be related to the rate of mean strain through an apparent 
turbulent or "eddy" viscosity given in incompressible flow by 
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Closure for the Reynolds heat flux-term, pCpTu] is treated in algebraic models by a form of Reynolds 
analogy. The Reynolds analogy is based on the similarity between the transport of heat and momentum 
and applied to the apparent turbulent conductivity in the assumed Boussinesq form 

pcpT'u'j = - 4 ) 

Experiments reveal that the ratio of the diffusivities for the turbulent transport of heat and momentum 
which is called the turbulent Prandtl number, Pr-f = prCp/kx. is a well-behaved function across the 
flow and in algebraic models is generally taken to be 0.9. Using the turbulent Prandtl number, the 
turbulent heat flux is related to the turbulent viscosity and mean flow variables as 


-pcpT'uj - 


Cppr dT 
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In the Baldwin-Lomax model the effects of turbulence are simulated by replacing the molecular 
coefficient of viscosity p.\i with the effective viscosity pm + pr in the stress terms of the laminar 
Navier Stokes equations. In heat flux terms k/Cp = p/Pr is replaced by pjPr + p T /Pr T - 
The Baldwin-Lomax model is a two-layer algebraic model in which p T is given by 


_ f ( pT)inner • V Vcrossover (Q _ g) 

l ( PT)outeri V — Vcrossover 

where y is the normal distance from the wall and y C rossover is the smallest value of y at wdiich values 
from the inner and outer formulas are equal. The eddy viscosity coefficient in the inner layer is based 
on the Prandtl mixing-length theory 

{pT) inner — pt ~ | a,’ | (C / ) 

The parameter t is the mixing length corrected with the Van Driest damping factor to account for the 
laminar sublayer 

£ — ky[l — exp(—y^ /A + )) (C — 8) 
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where k — 0.4. — 26. and \uj\ is the magnitude of the vorticity given by 
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The eddy viscosity coefficient in the outer layer is based on the distribution of vorticity which is used 
to determine the length scale and is given by 


{PT)outer = K Cep p F\ V AK E F K LEb(v) (C — 11) 

Fwake is found 

(C - 12) 

v-dif = (Vu 2 + ?; 2 + w 2 ) max — (\/ u~ + v 2 + w 2 ') min (C — 13) 

In equation C-13 the second term is taken to be zero (except in wakes). The quantities y ma x and F max 
are determined from the function 


where K = 0.0168 is the Clauser constant and Cep — 1.6 is an additional constant, 
via 
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where Cwk = 0.25 and 


F{y) = y | w |[1 - exp(—y^/A + ')\ 


(C - 14) 


In wakes the exponential term is set equal to zero. The quantity F max is the maximum value of F(y) 
that occurs in the profile and y max is the value of y at which F max occurs. The function F K leb is the 
Klebanoff intermit tency factor given by 


Fkleb{v) 



5.5 ( gj^gjj vr 1 

^ Vm ax ' . 


(C - 15) 


and Ckleb — 0-3. 

Although not seen explicitly, the Reynolds number enters into the computation of eddy viscosity 
through the computation of y~ ■ W hen the variables in equation C-10 are nondimensionalized. the 
following expression is obtained: 


y + 


Re p w r w _ 

\ r y 

V Mu' 


(C- 16) 


In high Reynolds number flows, as Re — ■ oo, the length scales l and F(y) in the inner and outer 
layers go to zero, which result in vanishing values of pr- For more discussion on the Baldwin-Lomax 
turbulence model, the reader is referred to the original paper (ref. 44). 
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